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Abstract 

F-theory GUT models favor a relatively narrow range of soft supersymmetry breaking 
parameters in the MSSM Lagrangian. This leads to the specific predictions that a 10 — 100 
MeV mass gravitino is the LSP, and the NLSP is quasi-stable, with a lifetime between a 
second to an hour. In a wide range of parameter space, the NLSP turns out to be a stau, 
though a bino-like lightest neutrahno is also possible. Focusing on F-theory GUTs with a 
stau NLSP, we study the discovery potential at the LHC for such scenarios. Models with 
a quasi-stable stau predict a striking signature of a heavy charged particle passing through 
the detector. As a function of the parameters of minimal F-theory GUTs, we study how 
many of such events to expect, and additional signatures correlated with the presence of 
quasi-stable staus. We also study the prospects for staus to become stopped in or near the 
detector, as well as potential ways to distinguish such models from minimal gauge mediation 
models with similar spectra. 
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1 Introduction 



A few decades of theoretical research has led to many novel ideas for what to expect at 
the TeV energy scale. Perhaps the most well motivated idea is that supersymmetry will 
soon be discovered. In addition to its aesthetic appeal, supersymmetry provides a relatively 
simple resolution of the hierarchy problem. There is also circumstantial evidence for su- 
persymmetry both from the prescient prediction of a heavy top mass, and the remarkable 
unification of the gauge coupling constants in the MSSM, which does not work with only 
Standard Model degrees of freedom. Though the lack of detection necessarily constrains the 
available parameter space, leading to a mild degree of fine-tuning (the mini-hierarchy prob- 
lem), supersymmetry remains among the most robust, and motivated options for beyond 
the Standard Model physics. 

Another motivation for supersymmetry is that at high energies it is predicted by string 
theory. Given the existence of meta-stable vacua where supersymmetry is broken, this 
suggests the attractive possibility that at least for matter charged under the Standard 
Model gauge group, supersymmetry persists down to the TeV scale. 

Even so, this by itself is not particularly predictive, because the details of supersymmetry 
breaking strongly affect the superpartner masses and interactions with the Standard Model. 
This in turn directly affects the expected collider signatures of a particular model. 

In the past few years, F-theory based models have emerged as promising candidates 
for particle phenomenology |T 48 . In these models, the Standard Model gauge group and 
matter originate in a compact region of the six extra dimensions of string theory. Physically, 
this corresponds to a regime where particle physics can be treated in a limit decoupled from 
gravity. Gauge groups localize on four real-dimensional subspaces, chiral matter localizes on 
two real- dimensional Riemann surfaces, and cubic interaction terms between chiral matter 
localize at points of the geometry. Flavor physics in particular appears to play a key role 



in constraining candidate geometries 14 , 24 , 32 , requiring a point of enhancement 32 



Combining this with the assumption that supersymmetry persists down to the TeV scale 
severely limits consistent scenarios. A deformation of minimal gauge mediated supersym- 
metry breaking seems to be the most natural option in F-theory GUTs [6||8||9], though 
other approaches are also possible]^ Parameterizing the effects of supersymmetry breaking 
through the vev of a GUT singlet (X) = x + in F-theory GUT deformations of min- 
imal gauge mediation models, the /i-term is generated through the Giudice-Masiero type 
operator 

^4^Xt_^^ (1) 



where A[/y is fixed by the internal scale of the geometry to be near the GUT scale. In 



^See for example 



49 



52 



and in the context of moduli mediation 4 53 



57 
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order to generate a weak scale value for /x, this requires a specific scale of supersymmetry 
breaking 

F ~ M^^akMcuT ~ 10^' GeV^ 

This fixes a narrow range of values for a/F ~ 10^ — 10^ GeV, which implies that the 
gravitino is the lightest superpartner (LSP) of mass 10 — 100 MeV, and that the next 
lightest superpartner (NLSP) is quasi-stable on timescales probed by collider detectors, with 
a lifetime on the order of one second to an hour. The primary decay chains of all MSSM 
superpartners involve the NLSP, which strongly influences the experimental signatures of 
an F-theory GUT. 

This begs the question: Which particle is the NLSP? In F-theory GUTs, constraints 
from other sectors of the model force the lightest stau ri to be the NLSP over much of 
parameter space. This is due to two effects: First, in minimal models with a point of 
i?8 unification, the messengers of the gauge mediation model tend to be in vector-like 
pairs transforming in the 10 © 10 of SU{5). In comparison with messengers in the 5 © 5, 
the higher dimension of the 10 representation lowers the mass of the stau relative to the 
gauginos. Requiring that gravity can be decoupled requires an asymptotically free SU{5) 
GUT group which constrains the number of vector-like pairs in the 10 © 10 to A'^io < 2. 
Moreover, there is also a stringy "PQ deformation" of the spectrum due to heavy gauge 
boson exchange of an anomalous U{1)pq |9|. This also has the effect of lowering the stau 
mass by an amount Apq which tends to be on the order of the weak scale. Even so, for a 
smaller range of parameter space a bino-like lightest neutralino Xi may be the NLSP. In a 
very particular Dirac neutrino scenario the right-handed sneutrino can also be the NLSP, 
but in this case the next to next to lightest superpartner (NNLSP) is also quasi-stable, and 
is typically the stau. 

If the NLSP is a stau, it will have dramatic consequences for the LHC. It is a heavy 
charged particle leaving the detector! Moreover, since R-parity is preserved, the number 
of such staus will always come in pairs. As a massive quasi-stable charged particle, the 
stau will initially register as a muon in the detector. The primary challenge is therefore to 
distinguish such "fake muons" from actual muons. Such stau events have extremely low 
Standard Model background, and can easily be isolated from actual muons through suitable 



selection cuts. Stau NLSP scenarios have been extensively studied in the hterature 58 63 



One of the most attractive phenomenological features of such scenarios is that because it 
is relatively easy to detect, it is possible to extract its mass. Correlating this with the 
four-momenta of other final state products, this provides a means to proceed up a decay 
chain, extracting the masses of other superpartners [62| . 

In this paper we perform a parameter space scan over F-theory GUTs, showing that 
over a broad range of parameter space the stau is an NLSP. Next, we study the prospects 
for detecting signatures of stau NLSP F-theory GUT scenarios at the LHC. See 22] for a 
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study of related F-theory GUT scenarios with a bino NLSP. The expected signatures in 
stau NLSP scenarios greatly depend on both the center-of-mass energy and the integrated 
luminosity, and with increases in each, the prospects for detecting staus improve. We find 
that with energies which will soon be achieved at the LHC, such signatures will indeed be 
within reach! 

Though the primary signature of stau NLSP scenarios is the presence of "fake muons" , 
the parameters of an F-theory GUT also strongly determine the additional signatures ex- 
pected. We find that the number of leptons in the final state also depends on the number of 
messengers A^io, and the size of Apq. For a single 10 © 10 pair of messengers and moderate 
PQ deformation, many of the staus are unaccompanied by other final states. For larger 
values of the PQ deformation, and for more messengers, additional sleptons become lighter 
than gauginos, providing new decay channels. This tends to increase the number of leptons 
expected with a given stau signature. 

We also study the prospects for generating staus which can become stuck in or near 
the detector and then decay at some later time. When these staus eventually decay to 
a gravitino and tau, observing such events would provide additional evidence for quasi- 
stable staus. Further, for staus stopped in the detector, correlating the presence of the tau 
produced from the decay of the stau with the initial collision event can in principle provide 
a means to measure the lifetime of the NLSP, and consequently the scale of supersymmetry 
breaking. 

Another remarkable feature of such scenarios is that because a stau NLSP is detectable, 
it is possible to reconstruct with high precision the masses of sparticles which decay into 
a stau. We review the expected level of precision in such mass fits, and then explain how 
this level of accuracy provides a means to break degeneracies between F-theory GUTs and 
minimal gauge mediation models with a stau NLSP which might otherwise be present with 
cruder signatures. 

The organization of the rest of this paper is as follows. In section |2] we review the 
salient features of F-theory GUTs. In section |3] we show that for a wide range of parameter 
space, the stau is the NLSP, while over a smaller range, the bino is the NLSP. We next 
move to collider signatures of F-theory GUTs with a quasi-stable stau NLSP. In section |4] 
we discuss details of the simulation performed and compute the production cross sections 
for supersymmetric events, and in section [5] we study search channels of potential interest. 
Section |6] discusses the prospects for observing stopped staus, and in section [7| we study 
whether mass reconstruction based on the stau NLSP can be used to distinguish F-theory 
GUTs from minimal gauge mediation models with similar spectra. Section [8] contains 
our conclusions. Appendix A contains additional details on the simulation. Appendix B 
contains some additional discussion on the selectron mass in F-theory GUTs and minimal 
gauge mediation, and Appendix C contains miscellaneous plots. 
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2 Parameter Space of F-theory GUTs 



Our aim in this paper will be to study the phenomenological consequences of a particularly 
minimal class of F-theory GUT models based on a point of Eg unification [32], and the 
corresponding supersymmetry breaking scenarios. See ji] for other supersymmetry breaking 
scenarios in F-theory GUTs, and their experimental signatures. 

In an F-theory GUT, all of the particle content of the Standard Model localizes on a 
spacetime filling seven-brane which wraps a complex two-dimensional surface. The matter 
fields of the MSSM localize on complex one- dimensional curves formed by pairwise inter- 
sections of the GUT stack with other seven-branes, and the Yukawa couplings of the MSSM 
localize at points of the geometry formed by triple intersections of seven-branes. In this 
type of setup, achieving realistic flavor hierarchies requires that these Yukawa points be 
close together. In minimal setups, there is a single point where the geometry enhances up 
to a point of Eg symmetry which is broken down to SU{6) and SO{10) along curves sup- 
porting respectively 5's and lO's of SU{5), and is broken down to SU{5)gut over the bulk 
of the seven-brane worldvolume. Minimal F-theory GUT models are especially predictive, 
and these are the scenarios we study in this paper. 

The supersymmetry breaking sector localizes on another seven-brane, which supports 
an anomalous U{1)pq gauge boson. This seven-brane intersects the GUT seven-brane, 
leading to the matter of the MSSM, all of which is charged under U{1)pq. Supersymmetry 
breaking is triggered by the vev of a GUT singlet X: 

{X)=x + e^F. (2) 

Supersymmetry breaking is communicated to the visible sector through gauge mediation 



64 -73 as well as through the exchange of a heavy anomalous U{1) gauge boson. See 74 
for a review of gauge mediation. In the deformation of minimal gauge mediation considered 
in [9] (see also Isl), the effects of supersymmetry breaking are transmitted to the messengers 



through the cubic interaction term: 

^ (i^d XYY\ (3) 



for messenger fields Y and Y' in vector-like pairs of complete SU{5) GUT multiplets. Once 
X develops a supersymmetry breaking vev, loops with gauge fields will then communicate 
this to the visible sector. Assuming that the interaction XYY' also localizes at the same 
point of Eg unification, it can be shown that in all but one Dirac neutrino scenario, the 



messenger fields are forced to transform in the 10 © 10 of SU{5)gut 32 . We denote 
the number of such vector-like pairs by Nio. Heuristically, the requirement that we can 
decouple gravity from the gauge theory requires that SU{5)gut is an asymptotically free 
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gauge theory, which imphes A'^io = 1 or 2. 

When the /z-term is generated through the higher-dimension operator of equation ([T]), 
the scale of supersymmetry breaking is fixed to be: 

F ~ M^eak ■ MgUT. (4) 

Introducing the characteristic scale: 

A = - (5) 

X 

of minimal gauge mediation, we can express the messenger scale Mmess in terms of x and 
the remaining parameters as: 

Mweak-McUT /„x 
X = Mmess ~ • (6) 

Phenomenologically, we require A ~ 100 TeV to achieve a realistic mass spectrum. Includ- 
ing all order one coefficients, we have [9]: 

F ~ 10^^ GeV^ (7) 
X ~ 10^2 GeV. (8) 

The specific range of supersymmetry breaking and messenger scales, as well as the con- 
strained messenger content constitutes one important prediction of this class of models. 
The vev x also doubles as the effective axion decay constant, which is an accord with the 
present window of viable values. 

At the messenger scale, localization on seven-branes allows an F-theory GUT to generate 
a weak scale value for the /^-parameter. Localization on seven-branes also implies that both 
B/j, and the A-terms vanish at the messenger scale: 

Bfl{Mmess) = (9) 
AijkiM^ess) = 0. (10) 

At lower energy scales the parameter Bfi is radiatively generated, and the ratio of the 



Higgs up and down vevs, Vu/vd = tan/3 ~ 20 — 30 is large 75 . Another consequence of 
generating S/i in this way is that proper electro- weak symmetry breaking then requires 
yuMi/2tan/3 > 0, where M1/2 denotes the gaugino masses (see for example 75 -77 ). Our 



sign convention is the same as 78 where M1/2, tan /3 > 0, which in turn implies > 0. To 
a large extent, the value of A determines the amount of fine-tuning in the Higgs sector. In 
particular, the current bound on the Higgs mass of ~ 114 GeV translates to a lower bound 
on A which we denote as Amm- The value of A^i^ also depends on the messenger content. 
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so that for A^'io = 1 and 2, we have: 



Amin(iVio = 1) = 50 TeV 
Amin(iVio = 2) = 28 TeV. 



(11) 
(12) 



To illustrate the broader features of this parameter space scan, we shall sometimes allow 
slightly lower values of A. The maximal value of A we consider is based on the requirement 
that we not overly fine-tune the Higgs sector, which to some extent is a matter of aesthetics. 
The upper bound on A we take is: 



A„ 
A. 



c(iVio 
c(iVio 



1) = 95 TeV 

2) = 53 TeV. 



(13) 
(14) 



Since the messenger scale also plays a role in setting the scale of the axion decay constant, 
we shall effectively treat Mmess as a fixed parameter. A change in A therefore corresponds 
to altering the scale of supersymmetry breaking. 

The soft masses at the messenger scale receive additional corrections from integrating 
out the heavy PQ gauge boson: 



L D — 47rapQexe<^ / d!^9 



(15) 



where ex and cj, respectively denote the PQ charges of X and the MSSM chiral superfield 
and M[/(i)pQ is the mass of the PQ gauge boson. Once X develops a supersymmetry 
breaking vev, this induces a shift in the soft masses at the messenger scale: 



m 



soft 



m: 



mGMSB 



2 

PQ- 



(16) 



Compatibility with minimal Majorana and Dirac neutrino scenarios 24 yields two dif- 
ferent choices for the PQ charges. In a normalization where the operator 



Anapqexex / d 9 



(17) 



U{1)PQ 

has the same coefficient in both Majorana and Dirac scenarios, the constant of proportion- 



ality q appearing in equation (16) is: 
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Figure 1: Spectrum of an F-theory GUT Majorana neutrino scenario with A'^io = 1, A = 50 
TeV for minimal (red left columns) and maximal (blue right columns) PQ deformation of 
order 200 GeV. At small values of the PQ deformation, the bino is the NLSP, whereas for 
larger values this transitions to the stau. Further note that at larger PQ deformation, the 
right-handed selectron and smuon are also lighter than the bino, and all of the sleptons are 
lighter than the second neutralino and lightest chargino. 



Technically speaking, there is slightly more freedom in the Dirac scenario, because there 
is also a U{1)b~l symmetry available in Dirac neutrino scenarios. Since the effects of the 
PQ deformation are similar in Majorana and Dirac scenarios, we shall focus on Majorana 
scenarios. 

The PQ deformation lowers the soft masses of the chiral matter fields, while raising 
the soft masses in the Higgs sector. For sufficiently large values of ApQ, this induces a 
tachyonic mode in the spectrum, which imposes an upper bound on the value of Apq on 
the order of 200 — 400 GeV, the precise value of which depends on A. Let us also note that 
there is also an experimental upper bound on the size of Apq because of the present 100 
GeV bound on the mass of the stau [79]. 

Summarizing, the parameter space of an F-theory GUT which we scan over is: 
• Fixed messenger scale Mmess ~ 10^^ GeV 
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• For A^io = 1 : 50 TeV < A(iVio = 1) < 95 TeV 

• For A^io = 2 : 28 TeV < A(A^io = 2) < 53.2 TeV 

• < ApQ < 200 - 400 GeV 

See figure [T] for a plot of tfie sparticle spectrum of a representative F-theory GUT witli 
A^io = 1, A = 50 TeV at minimal and maximal PQ deformation, and figure 23 of Appendix 
C for the case of A^io = 2 and A = 28 TeV. 



3 The LSP, NLSP and NNLSP 

The collider phenomenology of a supersymmetric model strongly depends on which sparticle 
is the lightest superpartner (LSP), the next lightest superpartner (NLSP), and sometimes 
the next to next lightest superpartner (NNLSP). In this section, we study the spectrum of 
F-theory GUTs, and show that the LSP is always the gravitino, and in Majorana neutrino 
scenarios the NLSP is either the bino or the stau. In Dirac neutrino scenarios, we find that 
for small values of the PQ deformation, a right-handed sneutrino is also a possible NLSP, 
though in such cases the NNLSP which is either a bino or stau is still quasi-stable on the 
same timescale as the NLSP. To perform our numerical analysis of the spectrum we use the 



program SOFTSUSY 78 



The narrow window for supersymmetry breaking implies that the gravitino is the LSP, 
with mass: 

fin F 

^3/2 = \hrTZ- ~ 10 - Mev. (19) 

V o Mpi 

The gravitino is stable, and all of the other MSSM superpartners eventually decay to it, 
plus their Standard Model counterparts. In this class of models, such gravitinos provide 



a viable dark matter candidate 17 , much as in the sweet spot model of supersymmetry 



breaking 80 



In many minimal gauge mediation models, the NLSP is a stau or a bino-like neutralino, 
and the situation in F-theory GUTs is no different. What is different in F-theory GUTs 
is that there is a specific range of values for the parameters of the model. Even without 
knowing the precise identity of the NLSP, the scale of supersymmetry breaking dictates 
that it will be quasi-stable, with a lifetime on the order of one second to an hour 12 , so 
that it will appear stable to a collider detector. 

Increasing the number of messengers or equivalently the dimension of the representation 
to which they belong tends to lower the masses of the scalars relative to the gauginos. For 
example, in a model with messengers in the 5 © 5, the scalar and gaugino masses at the 
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messenger scale are: 



'^scalar — V ^5 " C'scalar:;— A (20) 

<y 

'^gaugino -^5 ' C'gaugino A, (21) 

where schematically, the C's are constants fixed by the quadratic Casimirs of various rep- 
resentations, and the a's are the fine structure constants of the gauge groups under which 
the sparticles are charged. This is especially important for the present class of models, 
where the messengers are in vector-like pairs of 10 © lO's, the number of which we denote 
by A^io- The soft masses are proportional to the Dynkin index of the messengers. This 
means that each messenger pair in the 10 © 10 has the same effect as three in the 5 © 5, so 
that SA^io = A^5. 

The rest of this section is organized as follows. Besides the discrete parameter Niq, the 
low energy content of a minimal F-theory GUT is controlled by the messenger scale, the 
characteristic scale of minimal gauge mediation A and the size of the PQ deformation Apg. 
We first show that near the value of the messenger scale favored by F-theory GUTs, small 
changes in M^ess can switch the relative orderings of the bino and stau mass. After this, 
we study the effects of the PQ deformation. This latter deformation leads to qualitatively 
distinct mass spectra in the case of A^io = 1 models, but leads to more subtle changes in the 
relative masses when A'^io = 2. Next, we introduce some benchmark models, which crudely 
characterize these qualitatively different mass spectra. This is followed by a discussion of 
how shifts in the relative sparticle masses due to the PQ deformation affect decay topologies, 
which will be important later when we discuss candidate collider signals. 



3.1 Masses and the Messenger Scale 

We now discuss F-theory GUT scenarios at zero PQ deformation, but where we vary the 
messenger scale. For low scale gauge mediation models models with A^'io = 1, the stau is the 
NLSP. Increasing the messenger scale raises the mass of the scalars, because of the additional 
renormalization group evolution from the messenger to the weak scale. For example, in 
minimal gauge mediation models with A^io = 1, and A ~ 50 TeV, the bino becomes the 
NLSP for Mmess ~ 10^^ GeV. In minimal gauge mediation models with A'^io = 2, however, 
the stau remains the NLSP until Mmess ~ 10^^ GeV. 

Even when Niq = 1, the stau is the NLSP for a large region of F-theory GUT parameter 
space. Indeed, although the messenger scale tends to increase the mass of the stau relative 
to the bino, the PQ deformation lowers it back down. For example, in a Majorana neutrino 
scenario with A^'io = 1, Mmess ~ 10^^ GeV, A ~ 50 TeV, the bino and stau become equal in 
mass around Apg ~ 90 GeV. Figure |2] shows a plot detailing how Mmess and Apg determine 
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Logio(M„,ess/GeV) 




Figure 2: Contour plot indicating the relative masses of the bino, lightest stau and right- 
handed selectron as function of A pq and Mmess in the representative F-theory GUT scenario 
with A'^io = 1 and A = 50 TeV. There are three qualitatively different regions separated 
by the two lines where the bino is equal in mass to either the stau (upper line) or slectron 
(lower line). Here we also indicate the bino NLSP (upper left orange region) and stau NLSP 
(lower right cyan region) regimes of parameter space. 
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Mass (GeV) 




(GeV) 



Figure 3: Plot of the four lightest sparticle masses as a function of messenger scale M^ess- 
Here A'^io = 1, A = 50 TeV, and we have displayed the values for zero PQ deformation 
(upper lines) and Apg = 150 GeV (lower lines). Note that the bino mass remains constant, 
both as a function of messenger scale and of PQ deformation. 



the relative masses of the stau and bino. This also shows that for low enough messenger 
scales, and high enough PQ deformations, the right-handed selectron (and smuon, since 
they are degenerate in mass) can also become lighter than the bino. 

Figure [3] shows the masses of the four lightest Standard Model superpartners as a func- 
tion of Mmess and ApQ for a Majorana scenario F-theory GUT with A'^io = 1 and A = 50 



TeV. See figure [22] of Appendix C for the analogous plot with A^'io = 2, A = 28 TeV. Be- 
cause the effects of the messenger scale induce only logarithmic changes in the masses of 
the scalar partners, in the remainder of this paper we shall fix Mmess to be: 

Mmess = 10^' GeV. (22) 



3.2 PQ Deformation of the Masses 

In this subsection we study in greater detail the effects of the PQ deformation on the 
identity of the NLSP, and the other sparticle masses. As follows from the table of line (18), 



12 



the PQ deformation tends to lower the masses of all of the chiral matter superpartners 
relative to the gauginos. This effect is more pronounced for lower mass sparticles because 
the shift in the mass squared is: 



= m^o) + qA^Q (23) 



so that 



™ = ™(0)|l+^^ + o(^)). (24) 

Indeed, we find that the PQ deformation has a far greater effect on the masses of non-colored 
superpartners, so we focus on their masses in this subsection. 

Fixing the messenger scale, we now examine how the relative masses of the sleptons shift 
as a function of ApQ and A. Figure |4] shows the three qualitative crossing regions expected 
as one increases the value of the PQ deformation in F-theory GUTs with A'^io = 1 . The first 
crossover from a bino NLSP scenario to a stau NLSP scenario occurs at relatively low values 
of the PQ deformation. Moreover, by inspection of this plot, we observe that the value of 
ApQ for which this transition occurs is roughly constant as a function of A, going from a 
value of ApQ ~ 90 GeV for A ~ Amin to Apg ~ 100 GeV for A ~ A^ax- By inspection 
of the plot, there is a large range of parameter space where the stau is the NLSP, and a 
smaller sliver where the bino is the NLSP. Increasing Apq further causes the left-handed 
sleptons to become lighter than the second neutralino X2 lightest chargino Xi- This 
opens up new decay channels for these gauginos to left-handed sleptons, and significantly 
alters the decay topology. Further increasing Apg causes the right-handed selectron and 
smuon (which are mass degenerate) to become lighter than the bino. This alters the decay 
topology of these particles to three body decays mediated through an off-shell Xi- Figure 
5] shows a plot of the non-colored sparticles as a function of Apg for models with A^io = 1 
and A = 50 TeV. 

In A'^io = 2 scenarios, the stau is always the NLSP. Moreover, the right-handed selectron 
and smuon are also lighter than the bino, and all of the sleptons are lighter than the lightest 
chargino and second neutralino. As a consequence, there are no qualitatively different 
crossing regions for the relative masses as a function of Apq. On the other hand, the 
relative masses of the sparticles still depend on Apg, so that there will be some changes in 



the branching fractions. Figure 24 of Appendix C shows a plot of the non-colored sparticles 



as a function of Apg for models with A^io = 2 and A = 28 TeV. 
3.3 Benchmark Models 

In the previous sections we have studied how the mass spectrum depends on the parameters 
of an F-theory GUT. In the stau NLSP regime of A^io = 1 models, we have seen that there 
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A (TeV) 




Figure 4: As a function of Apg and A, the masses of the sleptons will change relative to 
the gauginos. For F-theory GUTs with A^io = 1 messengers, as Apg increases, there are 
three qualitative crossing regions where the stau becomes lighter than the bino, all sleptons 
become lighter than the second neutralino and lightest chargino, and the right-handed 
selectron and smuon become lighter than the bino. In the plot, the orange region to the 
left denotes the range of parameter space where the bino is the NLSP, and the cyan region 
to the right indicates the region of the stau NLSP. 
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Figure 5: Mass spectrum of the sleptons and lightest chargino and two hghtest neutrahnos 
with A^'io = 1, A = 50 TeV as a function of Apg. At values of Apg ~ 90 GeV, the 
stau becomes the NLSP. At somewhat larger values of the PQ deformation, the sneutrinos 
become lighter than X2 larger values, the right-handed selectron and smuon become 

lighter than the bino. The dashed grey line at 100 GeV denotes the present experimental 
bound on the mass of a quasi-stable stau. 
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are three qualitatively different orderings of the slepton masses relative to the electro-weak 
gauginos. For A'^io = 2 models, there is less variation in the mass spectrum, so that the 
primary qualitative differences are well represented by models with minimal and maximal 
PQ deformation. 

We now introduce representative "benchmark models" which reflect this behavior. Since 
the Dirac and Majorana neutrino scenarios are similar, we shall focus on the latter possi- 
bility. The benchmark models we consider are: 
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50 TeV, ApQ 


= GeV) 


(25) 
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= 175 GeV). 


(29) 



3.4 PQ Deformation of the Branching Fractions 



Shifting the relative masses of the sparticles also changes their branching fractions. In this 
section we study how the PQ deformation alters the branching fractions for particles which 
we expect to be produced by parton collisions. This includes gluinos and squarks, as well 
as second neutralinos and lightest charginos. 

Because the masses of the squarks are so high, the branching fractions for the squarks 
and gluinos are effectively constant as a function of Apq. By contrast, the PQ deformation 
has a significant effect on the non-colored sparticles, primarily because they are much 
lighter. To illustrate the effects of the PQ deformation, we first show the relative branching 
fractions of a few channels for Majorana neutrino scenarios for A'^io = 1 and 2, at minimal 
and nearly maximal PQ deformation for the benchmark models introduced in the previous 
subsection. With respect to these models, we compare the dominant decay channels and 
branching fractions of the gluino, lightest stop, lightest chargino and second neutralino. We 



use SDECAY 81 to compute the branching fractions numerically. 

First consider the decays of the gluino. Since the PQ deformation has only a small effect 
on the colored sparticles, there is only a few percent shift in the branching fractions: 



Channel 




Majlj\^ 


MajS 


Majjj'/ 


9 ^ tft^ 


40% 


33% 


22% 


21% 


g ^q^q^ 


27% 


33% 


36% 


37% 


g — b'fb^ 


21% 


21% 


16% 


16% 


g^b^b^ 


10% 


12% 


11% 


11% 



(30) 
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where here q refers to a first or second generation squark. Here we have omitted the decays 
to left-handed squarks, which only become relevant for the iVio = 2 case. Note that in going 
from A'^io = 1 to A'^io = 2 messengers, there is also a shift away from decays to stops and 
sbottoms in favor of first and second generation squarks. Nevertheless, the large branching 
fraction to stops in all cases suggest the presence of a large number of four top events. 



which should be relatively easy to see with low integrated luminosity 82 . 

Next consider the branching fractions of the stop. We similarly find that there is little 
shift for the two values of iVio and as a function of Apg: 



Channel 








M4nl 


f+ v° 

^1 ^ Al,2,3,4'' 


46% 


47% 


47% 


47% 




31% 


31% 


29% 


28% 


tt ^ xtb- 


22% 


23% 


25% 


25% 



(31) 



The branching fractions of the lightest chargino and second lightest neutralino exhibit 
stronger dependence on the PQ deformation. For example, the branching fractions of the 
lightest chargino are: 



Channel 
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11% 
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2% 


xt ''d^ 


0% 
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0% 
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0% 
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xt ^2^r 


0% 


7% 


7% 


10% 



(32) 



where here / denotes a first or second generation lepton. Similar A pq and A'^io dependence 
is also present in the branching fractions of the second neutralino: 
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0% 


8% 
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11% 



(33) 



Scanning over the A pq dependence of the branching fractions for the second neutralino 
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Figure 6: Dominant branching fractions of the hghtest chargino as a function of Apq with 
A^io = 1 and A = 50 TeV. As the mass of the sneutrinos becomes lower than the chargino, 
additional decay channels open up. 



and lightest chargino reveals that significant shifts occur for A^io = 1 F-theory GUTs when 
the sneutrinos become lighter than the chargino and second neutralino. Figures [6] and [7] 
show the respective changes in the branching fractions for the lightest chargino and second 
neutralino. Figures 25 and 26 in Appendix C show that for A'^io = 2, there is less of 
a dramatic shift in the branching fractions, in accord with the fact that all sleptons are 
already lighter than xf- 

Increasing Apg also increases the number of leptons present in the most likely decay 
chains. By inspection of all of the branching fractions, we see that as Apg increases, the 
first and second generation sleptons become light enough to become decay products for the 
lightest chargino and second neutralino. This decreases the branching fractions X2 ~^ t^t^ 
and xt ~^ '^i^T, meaning the expected decays will tend to have additional leptons. 



3.5 Right-Handed Sneutrino NLSP and Dirac Neutrinos 

In the above discussion we have assumed that the NLSP corresponds to a superpartner 
of a Standard Model particle. This is appropriate in Majorana neutrino scenarios, where 
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BRa2^X) 




Figure 7: Dominant branching fractions of the second neutrahno as a function of A pq with 
A'^io = 1 and A = 50 TeV. As the mass of the sneutrinos becomes lower than the second 
neutrahno, additional decay channels open up. 
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there are no additional light fermions beyond those already accounted for in the Standard 



Model. In F-theory GUT Dirac neutrino scenarios 24 , however, the Dirac mass is induced 
through the higher-dimension operator: 

^ ^ / d'eN^Nn (34) 

once Hd develops a supersymmetry breaking vev on the order of the weak scale. This type 
of F-term also induces an extremely small contribution to the mass of the right-handed 
sneutrino states, on the order of 0.05 eV. In certain ranges of parameter space, the vrS can 
correspond to the NLSP. In this case, the NNLSP is the stau or bino. Even so, this is rather 
inconsequential for collider phenomenology because of the highly suppressed interactions of 
the sneutrino with the MSSM. 

First consider the contributions to the mass of the sneutrino, all of which are correlated 
with supersymmetry breaking. Contributions from Planck suppressed operators will gener- 
ate a mass term on the order of the gravitino. Assuming this is the only contribution, this 
suggests the possibility that the vr could be the LSP, or at least very close in mass to the 
gravitino. 

Typically, however, the PQ deformation will generate larger mass terms. Indeed, since 
Nr has PQ charge q^Nr) = 3, there is a positive contribution of +3ApQ to the sneutrino 
mass squared. In the typical range of F-theory GUT parameter space where A pq is on the 
order of 100 GeV, ur is simply a heavy superpartner which is effectively decoupled from 
the MSSM. 

In principle, contributions from a U{1)b-l gauge boson can also contribute to the mass 
of the sneutrino, though this depends on assumptions about this sector of the model. For 
simplicity, we shall assume in this paper that the B — L gauge boson is heavier than the 
PQ gauge boson, but for the sake of argument, we now discuss what happens if we assume 
it interacts non-trivially with the supersymmetry breaking sector. Even when the B — L 
gauge boson has similar mass to the PQ gauge boson, it will have almost no effect on 



the sparticle spectrum. This is because the X field is not charged under U{1)b~l 32 
and so the analogue of the PQ deformation is necessarily suppressed by additional loops. 
When this gauge boson is lighter than the supersymmetry breaking scale, it can in principle 
induce additional contributions, much as in ordinary gauge mediation. In such scenarios, 
the mass of the sneutrino will then depend on the gauge coupling of U{1)b-l- Since all 
of the sleptons have the same magnitude U{1)b-l charge, the right-handed sneutrino will 
have mass lower than all of the other sleptons, and in particular the stau. Note, however, 
that the gauginos are not charged under U{1)b-Li so in principle either the bino or the 
right-handed sneutrino would in this case be the NLSP. 

Nevertheless, the effects of a right-handed sneutrino NLSP are fairly inconsequential for 
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the purposes of a collider study since the decay rate of other sparticles to the NLSP is at 
most comparable to the decay to the gravitino G. Indeed, the decay rate for a sparticle ip 
to decay to a gravitino plus its Standard Model counterpart ip is 

r. 

~ ■~- 771 

T{^l:^G + ^l:)r^J^, (35) 

with m the mass of tp. On the other hand, the decay channel provided by the right-handed 
sneutrino is suppressed by the scale Auv, and so assuming no additional loop suppression, 
will lead to a decay rate on the order of: 

T{^^un + .:)--j^. (36) 

In the context of an F-theory GUT where F/Auv is also on the order of the weak scale, 
it follows that decays to the right-handed sneutrino are at most comparable to decays to a 
gravitino. This implies that the NNLSP will still be quasi-stable. By abuse of terminology, 
we shall therefore refer to the NNLSP in such situations as the NLSP, since there is little 
distinction between these possibilities for collider studies. 



4 Event Generation of Long-Lived Staus 

In the previous sections we found that over much of the parameter space of minimal F-theory 
GUTs with Eg point unification, the stau is a quasi-stable particle on timescales probed by a 
collider detector. In the remainder of this paper we focus exclusively on Majorana neutrino 
F-theory GUT scenarios with a stau NLSP, and explore how associated stau signatures 
depend on the parameters of minimal F-theory GUTs. In addition to scanning over the 
parameters of F-theory GUTs, we also briefly comment on how the center-of-mass energy 
a/s affects the expected cross sections and distributions. Since the energy scales of operation 
7 TeV and 14 TeV are expected to be of particular relevance, we shall often focus on these 
cases. 

We have performed a Monte Carlo simulation of supersymmetric event production for F- 
theory GUTs at the LHC using PYTHIA [83], with basic detector effects included for various 
objects (see Appendix A for more details). The mass spectra and decay information are 



generated using SOFTSUSY 78 , SDECAY 81 and BRIDGE 84 , and this output is passed 



to PYTHIA in Susy Les Houches Accord format 85 . In subsection 4.1, we present the 
expected cross sections of F-theory GUTs as a function of the parameters A, Apq, as well 
as the center-of-mass energy ^/s. After the initial production of a supersymmetric event, 
the resulting sparticles will undergo a sequence of decays so that the final state contains 
an even number of staus (due to R-parity) which can be accompanied by Standard Model 
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particles. In subsection 4.2 we explain how we shall identify staus, taking into account 
some crude features of detectors. 



4.1 Production Cross Sections 



The number of expected supersymmetric events is determined by the production cross 
section of supersymmetric particles generated by parton collisions. Figure [8] shows the 
dominant leading order cross sections of an F-theory GUT with A^io = 1 and A ~ 50 TeV 



as a function of the center-of-mass energy ^/s of the LHC. See figure 28 of Appendix C 
for a similar plot for A'^io = 2 and A ~ 28 TeV F-theory GUTs. Figure |9] shows that the 
production of colored objects strongly increases by roughly two to three orders of magnitude 
in passing from ^/s ~ 5 TeV to y/s ~ 14 TeV. This is primarily because the mass of the 
gluinos and squarks are all above a TeV, so that the required energy transfer from partons 
must be at least of the same order of magnitude. Note, however, that even for smaller 
values of a/s, there is also significant production of non-colored sparticles such as lightest 
charginos, second neutrahnos, and Drell-Yan produced lightest staus. 
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Figure 8: Plot of the leading order cross sections of various production channels as functions 
of center-of-mass energy ^/s for the benchmark Maj^Jjj-, scenario. 

As a function of Apg, there is little change in the cross section of most channels. This 
is to be expected, because the percentage change in the mass of the squarks as a function 
of ApQ is small, and for the gauginos such as the gluino and lightest charginos and second 
neutralinos, the PQ deformation vanishes. The two notable exceptions to this are the Drell- 
Yan production of lightest staus and sleptons. This is because the mass of the lightest 
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Figure 9: Plot of the leading order cross sections of various production channels for F-theory 
GUTs with Nio = 1 as a function of A with a minimal value of ApQ such that the stau is 
the NLSP. Aside from Drell-Yan production of staus and sleptons, these cross sections are 



independent of Apg (see figure 10). 



sleptons more strongly depends on Apg, which in turn alters the expected production cross 
section. Figure [T0| illustrates that for an F-theory GUT with A^io = 1 and A = 50 TeV and 
as a function of Apg, Drell-Yan production of staus increases from 0.01 pb at moderate 
values of Apg ~ 100 GeV to roughly an order of magnitude more at Apg ~ 200 GeV, 
allowing it to become comparable to the production cross section of second neutrahnos and 
lightest charginos. Figure 29 of Appendix C shows a similar plot for the case A^io = 2 and 
A ~ 28 TeV. 

Next consider the dependence of the cross sections on the parameter A. In figure |9] we 
plot the leading order production cross sections in iVio = 1 models where a minimal value 
of Apg is chosen such that the stau is the NLSP. The production cross section for gluinos is 
initially larger than that of the non-colored sparticles, but then becomes smaller at larger 
values of A. See figure 29 of Appendix C for a similar plot of the case with A'^io = 2 and 
A = 28 TeV for vanishing Apg. A related issue is that as a function of A, the production 
cross sections have different slopes. This is because the parton distribution function (pdf) 
for the gluons and sea quarks are different as a function of the momentum transfer fraction 
X. In particular, it is known that the gluon pdf significantly decreases at large x. Thus, as 
the sparticles increase in mass, the resulting momentum transfer from partons must be that 
much larger, leading to a pronounced decrease in sparticles produced from gluon colhsions 
rather than quark collisions. At a theoretical level, increasing A also tends to lead to more 
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Figure 10: Plot of the leading order cross sections of the dominant production channels 
functions of Apg for A^'io = 1 and A = 50 TeV. 
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fine-tuning in the Higgs sector. In other words, the more motivated possibihty of lower A 
is also the one most accessible to the LHC, and this is the case we shall typically focus on. 



4.2 Candidate Staus 

After the production of a supersymmetric particle, it will undergo a sequence of cascade 
decays down to the quasi-stable stau NLSP. As opposed to many supersymmetric scenarios 
with either a quasi-stable bino or wino, which register as missing transverse energy, a stau 
will leave a more prominent signal, because it is electromagnetically charged. 

Though a complete detector simulation is beyond the scope of the present paper, the 
detector resolutions and dp for the velocity (3 = v/c and momentum p of candidate staus 
are ^Oj (see also |86]): 

(7/3 ^ (T„ kip , / m? k-i GeV , , 

^ = 0.028/3, -J^ = -^(Bk2\ 1 + ^®- 37 

p p GeV Up p 

where ki = 0.0118%, ^2 = 2% and k^ = 89%. Since staus are quasi-stable, and have the 
same charge as muons, many will pass the level one (LI) and level 2 (L2) muon triggers, and 
will initially register as muons, provided they have high enough velocity. Indeed, standard 
data collection and reconstruction of charged tracks in the inner tracking chamber and muon 
detection system will record many muon events, of which some will actually correspond to 
staus. 

The staus generated by cascade decays from heavier particles will typically have large 
velocities. On the other hand, nearly all muons will have large /3. The Standard Model 
background for stau candidates mainly comes from muons whose velocity /3 is inaccurately 
measured. The typical Standard Model processes involving muons are dominantly the 
single boson production W/Z + jets, as well as top and bottom pair production ti and bb. 
With the cut (3 < 0.91 on stau candidates, a rejection rate of ~ 1000 can be achieved on 
background muonsj^ In order to differentiate muons from staus, we shall therefore impose 
the selection cuts for candidate staus: 

• 0.67 < 13 < 0.91 

• Pt> 20 GeV, \ri\ < 2.5. 

Let us comment on the selection cut for (3. As already mentioned, the upper bound on 
P is necessary to reduce contamination from actual muons. The lower bound reflects the 
requirement that candidate staus which move too slowly will not be properly identified 



Using fast-moving staus is also possible, see 87 for a recent study. 
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with the correct bunch crossing. At the ATLAS detector, for example, the muon trigger 
efficiency drops rapidly for slow moving particles with velocity f3 < 0.8 in the barrel, and 
(3 < 0.7 in the endcap of the detector 88 . In principle, the bound (3 > 0.67 can be relaxed 
to 0.6 |88|, though the detector efficiency is very low. Lowering the cut in this way requires 
that slow-moving staus can be reconstructed through an offline analysis using the monitored 
drift tube (MDT) data. In section [s] we shall discuss how inclusive stau signals depend on 
this lower bound. For /3 > 0.8, the efficiency is very close to 100%. The cuts on pt and t] 
are included because this is the acceptance range for the detector. 

Figure [n] shows a plot of the (3 distribution of staus for the benchmark Maj[J|j~, scenario 
for 10, 000 supersymmetric events at ^/s = 7 TeV and 14 TeV with the selection cuts on 
and \ri\ imposed. This plot illustrates that the (3 distribution is strongly peaked near /3 ~ 1, 
and is much lower at smaller values of (3 near the lower bound of 0.67. Another interesting 
feature of this plot is that increasing the center-of-mass energy from ^/s = 7 TeV to 14 
TeV tends to push up the number of high /3 staus. This means that once the selection cut 
on P is imposed, the acceptance of staus will drop slightly for the 14 TeV case. Of course, 
this is counterbalanced by the fact that the overall number of supersymmetric events is 
expected to be higher at y/s = 14 TeV. Said differently, because the distributions have 
roughly the same shape for 10, 000 events each, the much higher number of supersymmetric 
events expected at ^/s = 14 TeV will increase the signal, and the discovery potential. 
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Figure 11: P distribution of candidate staus generated with 10,000 supersymmetric events 
for the benchmark Maj^jj-, scenario at center-of-mass energies ^/s = 7 TeV and 14 TeV. 



Not all staus, however, will pass the muon trigger system. For slow enough staus, the 
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timing of a "muon event" may not correlate with the bunch crossing of the colhder beam, 
and will therefore not pass the initial trigger system, and must instead be detected through 
a dedicated offline analysis. Though the direct detection of slow-moving staus is therefore 
more challenging, their low velocity will allow some of them to become stopped in or near 
the detector, providing a means to measure additional properties of the stau. We shall 
return to the prospects for such stopped stau scenarios later in section |6} For this reason, 
in the next section we focus on scenarios where candidate staus pass the muon trigger. 
Indeed, because such staus are produced through the decay of heavier sparticles, a large 
fraction of them will have sufficiently high velocity to pass the LI and L2 muon triggers. 

5 Search Channels 

We now discuss potential search channels for stau NLSP scenarios in F-theory GUTs. In 
many supersymmetric extensions of the Standard Model, a common signature consists of 
searching for cascade decays which terminate with a large missing transverse energy carried 
by either the LSP or quasi-stable NLSP, accompanied by some number of Standard Model 
states. Here, this extra energy carried by the analogous quasi-stable supersymmetric state 
is more visible, since all cascade decays terminate with the stau NLSP. 

The characteristics of a long-lived charged particle passing through the detector would 
already provide a striking deviation from the Standard Model. Combined with other evi- 
dence, such as a determination of the spin of the particle, this would provide evidence in 
favor of a stau NLSP scenario. It is therefore important to determine how many such stau 
events to expect at the LHC, and how the number of events depends on the parameters of 
F-theory GUTs. The background from Standard Model processes is quite low in such sce- 
narios, and is mainly due to improperly measuring the velocity of muons. We shall return 
to specific selection cuts which can be used in various search channels, and how these can 
be used to further reduce Standard Model (SM) background. 

One of the most promising discovery channels for staus is based on summing over all 
possible candidate events in inclusive channels which contain one or two stausj^ We study 
how the parameters A and Apg of an F-theory GUT, as well as the parameters of the LHC 
such as the integrated luminosity and center-of-mass energy affect the discovery potential 
of finding staus at the LHC. We find that though the total cross section strongly depends 
on A, there is little dependence on Apg. This means in particular that such search channels 
can be used as a means to constrain the value of A. After this has been done, other search 
channels can then be used as a means to constrain Apg. 

More refined information about the parameter dependence of F-theory GUTs can be 

^The reason we also consider one stau events is that the other stau expected based on R-parity may not 
pass the selection cut imposed to differentiate staus from muons. 
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extracted from exclusive channels which require a specific combination of particles beyond 
the ubiquitous staus. Indeed, one particularly clean signature consists of just one or two 
staus, and nothing else. Such signatures can arise from Drell-Yan pair production of lightest 
staus, or alternatively lightest charginos. In the latter case, some missing energy will also 
be carried away by neutrinos. Since such signatures are most directly sensitive to the mass 
of the stau, they provide an especially clean means to track the dependence on Apg. Other 
exclusive channels with low background are based on stau events with some number of 
additional leptons. 

Even when we cannot distinguish the stau from a muon, the corresponding events with 
additional leptons are difficult to mimic with Standard Model processes. Such events cor- 
respond to leptons with the same charge which we refer to as same-sign leptons, as well as 
events with three or more leptons. 



5.1 Inclusive Stau Channels 

In this subsection we discuss inclusive stau channels obtained by summing over all signals 
containing one or two staus. The inclusive > 1 stau signal has the largest cross section and 
therefore provides the best channel for a discovery. Even though all supersymmetric events 
will contain two staus because of R-parity, imposing selection cuts used to isolate a signal 
will lead to some loss of signal. In particular, some two stau events may instead register 
as one stau events, or may not register at all. Letting (jf^^^ denote the total production 
cross section of supersymmetric events, we define the acceptance ei,2f of one and two stau 
events as the ratio of the observed cross section of one and two stau events a^l'^ and cr"!** to 

^SUSY. 
^tot 

ohs 

— ^SUSY V'^°/ 

for i = 1,2. In the remainder of this subsection, we analyze how sensitive the observed cross 
sections are to the choice of (3 cut imposed, and then discuss the sensitivity of this cross 
section to the parameters A and Apg. After this, we discuss the prospects for discovery in 
this channel after further selection cuts are imposed. 



5.1.1 /3min and Cross Sections 

As mentioned in section |4| the number of expected staus depends on the choice of selection 
cuts, and in particular, the lower bound f3min on the velocity (3. To a certain extent, this 
lower bound depends on how well slow-moving staus can be identified by an offline analysis. 
To gauge how much loss of signal to expect, we compute how (T°|'' and ctjI* depend on (3min- 



In figure 12 we plot this dependence for the benchmark Majillj-, model, where we also show 



the acceptances eif and ■ This plot illustrates that for ^/s = 14 TeV, the cross section 
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of one stau events a°'^^ is always a few times larger than that for two stau events o'2f'^. 
This is primarily because a large portion of staus have velocity /3 > 0.91, and are therefore 
eliminated by the selection cut. 



Figure 12 also illustrates that there is little dropoff in the value of a^^^ or ctj^* as a 
function of f3min until f3min becomes comparable to 0.7, which is around the value of the 
selection cut taken anyway. One curious feature of this plot is that in the case where ^/s = 7 
TeV, the one stau cross section cr"^"* is almost the same as a2f^ when f3min is small, but then 
increases as f3min increases up to ~ 0.7. This occurs because as the cut becomes more 
stringent, some of the previously accepted two stau events will now be rejected, and will 
instead be recognized as single stau events. Comparing the total acceptance at y/s = 7 TeV 
and 14 TeV, note that although the total cross section at 7 TeV is lower, the acceptance is 
higher. This is again because the stau velocity tends to be smaller at the lower center-of- 
mass energy, meaning that fewer events are rejected based on the selection cut (3 < 0.91. 



5.1.2 A and Apg Dependence 

Having seen that in the case of representative benchmark models that there is little change 
in the observed cross section as we change the lower bound on the selection cut (3min < /3, 
we now fix the value of this cut at f3min = 0.67, as in section |4} In this subsection we 
focus on how the production cross sections depend on the F-theory GUT parameters Niq, 
A and Apg. Figure |9] shows that as we increase the parameter A, we find that the total 
production cross section atot decreases. This is to be expected because increasing A also 
increases the mass of all of the sparticles. Indeed, even though increasing Apq will tend 
to lower the mass of the sleptons and in particular the stau, note that in such a situation, 
the mass splitting between the heavier sparticles and stau will become more pronounced. 
This in turn means that staus generated by a cascade decay will tend to be more energetic, 
having higher velocity, and will therefore not pass the selection cut on /3 as easily. 



Figure [13] shows that at least for these inclusive stau signals, the dependence on the 
PQ deformation Apg is quite mild. This is mainly because the bulk of produced staus are 
generated from decays of colored sparticles which have masses largely independent of this 
deformation. If stau signals are observed, this inclusive channel can then be used as an 
effective means to constrain the parameter A, once a fit to the F-theory GUT parameters 
has been performed. 



5.1.3 Discovering Staus 

Imposing additional selection cuts on inclusive stau channels can further eliminate contam- 
ination from background muons produced through Standard Model processes, leading to 
a cleaner signal and a better chance of discovery. In addition to the basic cut (3 < 0.91 
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Figure 12: Observed cross sections and acceptances for events with one and two stau 
candidate events as functions of the lower velocity cut /3mm for the benchmark model 



defined in subsection 



3.3 



The acceptance is by definition proportional to the 



observed cross section (see equation (38)). Also displayed is the total number of one and 
two stau candidate events. The upper and lower panels are for y/s = 7 TeV and 14 TeV 
respectively. 
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Figure 13: Plots of the observed cross sections of inclusive one and two stau signatures at 
a/s = 14 TeV as a function of ApQ for stau NLSP F-theory GUT scenarioss with A^'io = 1 
and A = 50 TeV. 



mentioned in subsection 4.2 which we use to avoid contamination from muons, we shall also 



include the event selection cuts in 60,61 :^ 

• At least one hadronic jet with px > 50 GeV and a calorimetric EJ^^^'^ > 50 GeV 
(trigger requirement); 

• me// > 800 GeV, 

where fUeff is the total invariant mass of the event constructed starting from the transverse 
momentum of the high px jets and muons (or muon-like particles) 



m,in{A,Njst) 



1=1 



Pt + 



E 

i=l 



Pt 



(39) 



The Standard Model background after all these cuts is estimated to be affP'^ ^ 1 fb and 



has only a mild effect on the signal acceptance 60,61 . Note that since this cut requires 



Though we do not include it here, additional discrimination from background muons can be achieved 



as m 



60 



61 



by requiring p > rrif-^ 



/3-0.05 



where mf^ is obtained from the stau mass measurement 



(see section 7^. This cut significantly increases the level of muon rejection, but only affects candidate stau 
events at the level of a few percent. 



31 



at least one hadronic jet, it will reject events generated by Drell-Yan production of staus 
and charginos. In practice, this means that such signals can be treated as independent 
constraints on the parameter space of F-theory GUTs. 




40 50 60 70 80 90 

A (TeV) 

Figure 14: Contour plot of the inclusive > Ifi signal in parameter space A and ApQ for 
^/s — 14 TeV and Aio = 1 messengers. The three contours correspond to the value of the 
cross section at 0.5 pb, 0.05 pb, and 0.01 pb. To detect five inclusive one or two stau events, 
this amounts to an integrated luminosity of respectively 10 pb~^, 100 pb~^ and 0.5 fb~^. 

The estimated upper bound of (t'sm'^ i$ 1 means that even for 1 fb~^ of integrated 
luminosity, we can expect around one Standard Model event. Assuming one Standard 
Model event as background, a 5(7 level discovery would then correspond to roughly seeing 
five times as many candidate stau events. To gauge the necessary integrated luminosity 
required to claim discovery, we have computed the inclusive cross section at ^/s = 14 TeV 
of one and two stau candidate events as a function of Aio, A, and Apg by imposing the 
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further selection cut mentioned above. In figure 14 we plot for A'^io = 1 models the resulting 
contours for cross sections of 0.5 pb, 0.05 pb, and 0.01 pb. See figure 27 of Appendix C for 
the analogous plot for A^'io = 2 models. To detect five inclusive one or two stau events, this 
amounts to an integrated luminosity of respectively 10 pb^^, 100 pb~^ and 0.5 fb^^. Note 
that this in some sense overestimates the number of Standard Model background events, 
which upon allowing for fractions of an event for each integrated luminosity is roughly 0.01, 
0.1 and 0.5. 

At lower center-of-mass energies ^/s ~ 7 TeV, the total production cross section is 
smaller, so the prospects for discovery are lower. Even so, for large enough integrated 
luminosities, it may still be possible to discover such stau events. Indeed, returning to 
figures |8] and 12, we see that the total cross section of inclusive stau events at 7 TeV is 
lower by a factor of 10 compared with the cross section at 14 TeV, so that the required 
integrated luminosity for discovery is roughly a factor of 10 higher compared with what 



is required in the 14 TeV plots of figure [14] and figure [27] of Appendix C. For example 
with a few hundreds of pb~^ data delivered at 7 TeV, we can expect to discover staus with 
A = 50 TeV for A^io = 1 or A = 28 TeV for A^io = 2. Depending on how long the LHC 
operates at 7 TeV, higher values of A may also be within reach. 



5.2 Exclusive Two Stau Channel 

Additional information about the parameter dependence of F-theory GUTs can also be 



we 



extracted by examining more specialized channels. In particular, in subsection 5.1.3 
imposed a selection cut on inclusive stau signatures requiring at least one hadronic jet. 
This excludes, however, signatures such as Drell-Yan production of staus and charginos, 
which can lead to two opposite sign staus. In the case of Drell-Yan produced charginos, the 
decay of the charginos will also generate neutrinos which will leave the detector as missing 
energy. Because the stau mass depends strongly on Apg, this correlates the expected cross 
section with ApQ. Moreover, with the additional information on the measured stau mass 
(see section [T]), the stau production cross section can also be used to determine the spin of 



the stau 89 . Though the specifics of the event selection are different, this same channel 



has been used in Tevatron searches for long-lived charged particles 90,91 . 

In this subsection we consider events with only two staus, and no additional (detectable) 
Standard Model states. More explicitly, we require there is at least one stau candidate which 



passes the cut on /3 and rj introduced in subsection |4.2[ and which passes the slightly more 
stringent cut that pt > 100 GeV for any stau candidate. For candidate events in which 
only one stau passes this cut, we shall still require that an additional "lepton" has px > 100 
GeV. To exclude the presence of final states with additional particles, we also veto events 
with jets satisfying pt > 20 GeV and I77I < 3 and extra leptons satisfying pt > 5 GeV 
and |?7| < 2.5. To further distinguish between staus generated by Drell-Yan production of 
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staus and charginos, it is also helpful to impose the cut i?™'*'^ < 10 GeV which eliminates 
staus created from charginos which decay to a highly energetic neutrino plus a stau. In the 
following we consider exclusive two stau events with and without this cut. 

We now study the dependence of the resulting exclusive two stau events as a function 
of ApQ, fixing a particular representative value of A = 50 TeV for A'^io = 1 and A = 28 TeV 
for A'^io = 2 F-theory GUTs. In figure 15 we plot this dependence. First consider the case 
of A^io = 1 F-theory GUTs. When the further cut EJp^^^ < 10 GeV is not imposed, there is 
a two-fold degeneracy in the expected cross section. This is because the two contributions 
from Drell-Yan produced staus and charginos behave differently under shifts of the PQ 
deformation. Whereas the PQ deformation decreases the mass of the stau, and thus leads 
to a monotonically increasing contribution from Drell-Yan produced staus, for Drell-Yan 
produced charginos, the overall production cross section of charginos is independent of 
ApQ, but as seen in section |3| the subsequent branching fraction of charginos to a stau and 
neutrino decreases at larger ApQ, as new decay channels become available. This has the 
effect of lowering the contribution to this channel. In tandem with the contribution from 
Drell-Yan produced staus, the net effect is to produce a small dip in the cross section as a 
function of Apg. Note, however, that once the further cut i?™'^'* < 10 GeV is included, the 
cross section is again a monotonically increasing function of Apg, as expected. 

Next consider the case of A'^io = 2 F-theory GUT models. In section [3] we found that 
the change in the chargino branching fraction is much milder in such scenarios, essentially 
because no new decay channels open up as Apg is changed. As a consequence, even without 
imposing a further EJp^'^'^ < 10 GeV cut, we see that the exclusive two stau channel is a 
monotonically increasing function of Apg. Imposing this cut leads to a shift of about 0.005 
pb in the overall cross section, which is largely independent of Apg. 

Since the stau is the sparticle with mass most sensitive to the PQ deformation, it 
is especially interesting to use the expected number of events to extract more detailed 
information on the PQ dependence. Imposing the further selection cut EJp^'^^ < 10 GeV, we 
note that in figure 15, the overall cross section ranges from roughly 0.005 pb at low Apg 
in the Niq = 2 scenario, to about 0.02 pb in both A^io = 1 and A'^io = 2 scenarios. For 10 
fb~^ of integrated luminosity, this translates to 50 — 200 expected events. 
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Figure 15: Cross sections of exclusive events with only two staus as functions of ApQ in 
stau NLSP F-theory GUTs. The upper panel is for N^^ = 1 and A = 50 TeV, while the 
lower one is for iVio = 2 and A = 28 TeV. The upper dashed line corresponds to the signal 
without a cut on E^jl'***, while the lower solid line corresponds to the signal with < 10 

GeV. 
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5.3 Two Stau Plus Lepton(s) Channels 



In this subsection we study two stau events which include one or two leptons. Such final 
states typically originate from the decay of the second neutralino and lightest chargino as 
in the decay processes: 

Xl ^ P + ^fi^ + r^ + /+ + /-, 

f± + 

xt t'i + ^T, 

Dr + T^^ + + r+ + 

where as usual in collider studies, we have distinguished the somewhat different r and other 
leptonic final states (Z's). r's which decay leptonically to e and will contribute isolated 
leptons to the final state. Given the fact that the lepton efficiency is high, the relative size 
of different signatures with different number of leptons encodes the information about the 
branching factions of chargino and neutralino decays, and therefore can be used to probe 
the size of ApQ, as well as to serve in reconstructing the mass of the second neutralino and 
lightest chargino, a point we shall return to in section [7| 

Here we consider the following signatures: 

• 2 fi + (> 1/) + jet 

• 2 fi + (> 1/) 

• 2 fi + (> 21) 

where the 2f events are selected as in the previous subsection, and "0 jet" refers to events 
which are vetoed if they contain a jet with pt > 20 GeV and \r)\ < 3. For these channels, 
the background is again suppressed by the strong stau selection cuts, as well as by requiring 



additional leptons. The dependence of these signals on Apg is shown in figure 16 for F- 
theory GUTs with A^io = 1 and A = 50 TeV, as well as A^io = 2 and A = 28 TeV. By 
inspection, we see that the inclusive stau signals with leptons and jets have large cross 
section as expected and increase monotonically with Apg. 

By contrast, the stau signal with leptons but no jet has relatively small cross section. 
This is also expected because such signals are produced entirely by electro-weak processes, 
and so are not connected with the decay of a squark or gluino, which have much higher 
production cross sections when compared with charginos and neutralinos. In addition, note 
that this signal has very little dependence on Apg, which is to be expected because the 
gaugino masses are independent of Apg. 
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Figure 16: Plots of the cross sections of candidate events with two staus and > 1 lepton + 
jets, > 1 lepton and > 2 leptons as functions of Apg in the stau NLSP regime of F-theory 
GUTs with iVio = 1 and A = 50 TeV (upper panel) and A^io = 2 and A = 28 TeV (lower 
panel) . 
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5.4 Inclusive Lepton Signals 



In the above analysis we have focussed on signatures which aim to distinguish staus from 
muons. Even when such staus cannot be distinguished from muons, they can still generate 
signatures with low Standard Model background. Since Standard Model processes copiously 
produce leptons, such signatures require either the presence of multiple leptons, or specific 
charge combinations of leptons combined with a large px cut. 

In this subsection we study inclusive lepton signals where wc do not explicitly attempt 
to differentiate staus from muons. In other words, we now consider any stau as a "muon" 
if it satisfies: 

• /3> 0.67 

• Pt> 20 GeV, \ri\ < 2.5 

With respect to these cuts which allow fast-moving staus, the inclusive lepton signatures 
we consider are: 

• SS: A pair of same-sign isolated leptons. 

• 3/: Three isolated lepton candidates. 

• 4:1+: Four or more isolated lepton candidates. 

In order to reduce background, we also require these signal events to have: 

• At least two hard leptons with pr > 100 GeV 

• At least two hard jets with px > 150 GeV. 

Since staus passing the "muon" cut should always have pr ^ 100 GeV, the px cut on the 
two hardest leptons only mildly affect the signals, but can significantly reduce Standard 
Model background. For the SS signal, the Standard Model background is negligible. For the 
3/ signal, the primary source of Standard Model background is from WZ and tt production. 
We reduce the WZ background by requiring that no pair of charged leptons reconstructs to 
an invariant mass mz ± 10 GeV. The tt background is potentially large if the muon from a 
b quark decay is isolated. To suppress it, we require two hard leptons with pt > 100 GeV 
to reduce the acceptance of leptons from W decay, and additionally require two hard jets 
with Pt > 150 GeV. Alternatively, if one is interested in the signal from electro-weakino or 
slepton productions, one can also use a jet veto to suppress the ti background, which we 
will not discuss here. For the 41+ signal, the Standard Model background is pre-dominantly 
from ZZ production. We reduce this background using the same type of cut as in the 3/ 
signal. 
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Figure 17: Plots of the cross sections of inclusive lepton events in which fast-moving staus 
are identified as leptons for two F-theory GUT scenarios with A^io = 1 at A = 50 TeV 
(upper panel) and A = 75 TeV (lower panel). The same-sign lepton (SS) tri-lepton (3/) 
and at least four lepton (4/-I-) signals respectively correspond to roughly increasing and 
decreasing functions of Apg. 
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Figure 18: Plots of the cross sections of inclusive lepton events in which fast-moving staus 
are identified as leptons for two F-theory GUT scenarios with A'^io = 2 at A = 28 TeV 
(upper panel) and A = 42 TeV (lower panel). The same-sign lepton (SS) tri-lepton (3/) and 
at least four lepton (4/-I-) signals respectively correspond to roughly increasing, constant, 
and decreasing functions of Apg. 
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Let us now discuss the primary origin of these signatures and our theoretical expecta- 
tions for the corresponding dependence on the PQ deformation. For the same-sign (SS) 
isolated lepton signal, gluinos and squarks can decay to pairs of same sign lightest charginos 
xfxf pairs of second lightest neutralinos X2X2- The subsequent decay of the light- 
est chargino for example via xf ~^ '^t + and the second neutralino for example via 
X2 — 7- Tf*" + with the r decaying hadronically, will generate SS lepton signals. 



As can be seen from figures 17 and 18, the SS signal is roughly a decreasing function of 
ApQ. To see how this comes about, note that in the case of same sign staus generated by 
lightest chargino and second neutralino decays, as A pq increases, the branching fraction to 
states without additional charged leptons decreases. Note that increasing the value of A pq 
simply diverts these decays through channels which generate additional charged leptons in 
the final state. 

Tri-lepton (3/) and four or more lepton (4/+) signatures correspond to processes involv- 
ing two fast-moving staus, as well as additional e or /i leptons, and can be generated in a 
variety of different cascade decays. For example, focussing on the contribution from a single 
branch of a cascade decay, the 3/ signal can be generated by decays involving — )■ / + z/, 
Xi — X2 ~^ + or other channels where X2 decays to a stau and some taus. Note 
that we can get contributions to the 3/ and 41+ signatures depending on whether the taus 
in the final state decay hadronically or leptonically. Figure 17 shows that for A^io = 1, both 
of these signals are roughly increasing functions of ApQ. As Apg increases, the masses of 
gl, a*!,, T2, Ue-i ^fj. and Uj- decrease and eventually drop below the mass of X2 and xf^ while 
en and Jin both become lighter than x'j'j^ Consequently, the branching ratios of X2 and xf 
into these states increases. Since these states typically decay to ri with additional leptons, 
the increase of Apg also increases the 3/ and 41+ signals. 



In the case of A'^io = 2 models, figure [18] shows that the 31 signal is roughly constant, 
while the increase in 41+ signatures becomes more prominent than for A'^io = 1. The 
difference in the A^io = 1 and A'^io = 2 cases is due to a subtle tradeoff in the increase and 
decrease of the branching fractions of heavier sparticle states. For example, the branching 
fraction of x^'s via — )■ / + increases as Apg increases, while the branching fraction of 
X2 + decreases. For iVio = 1, the net contribution from these effects leads to an 

increase in the overall number of 3/ signals, while for A'^io = 2, the positive and negative 
contributions to this branching fraction more exactly balance out, leading to a constant 
ApQ dependence. 



^Note that in the case of A^io — 2 F-thcory GUTs, this statement already holds for vanishing PQ 
deformation. 
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6 Stopped Staus 



Up to this point, we have treated the stau NLSP as stable on timescales probed by coUiders. 
In this section we study the prospects for F-theory GUT stau NLSPs to become stopped in 
or near the detector, and to decay at some later time to a tau and a gravitino. Searching 
for stopped staus through their decay products would provide additional evidence that the 
quasi-stable massive charged particle is in fact a stau. When the origin of the stopped stau 
can be reconstructed, observing its decay would provide a measurement of its lifetime, and 



consequently the scale of supersymmetry breaking 86 , 92 95 



An analysis including all details of the ATLAS and CMS detectors, and in particular the 
ways that they are situated in their respective caverns is beyond the scope of the present 
paper. To extract the main physical features of stopped stau scenarios, we consider an 
idealization where after the initial production of staus, these particles enter a large spherical 
ball of either iron or carbon. Using the overall production cross section of staus as well 
as the associated velocity distribution obtained in section |4| we simulate 100 fb~^ of LHC 
data and compute the total number of stopped staus N{R) contained within a distance R 
from the collision point. Here we focus on candidate staus with px > 20 GeV and \t]\ < 2.5. 
As expected, we see that integrating over an ever larger ball leads to more stopped staus. 
Though small, the number of such stopped staus is not completely negligible. In figure 
we show a plot of N[R) for stopping in iron and carbon for the Maj^^j^, scenario. In 
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this plot, we see that the expected number of stopped stau events is on the order of a 
few thousand after passing through ~ 5 meters of either iron or carbon. More generally, 
increasing A tends to decrease the total production cross section of staus, and thus the total 
number of staus which could become stopped. 

The effects of the PQ deformation on the number of stopped staus is more subtle. Fixing 
the mass of the other sparticles while lowering the mass of the stau will cause the staus 
which are produced through a decay chain to be more energetic, which in turn skews the 
velocity distribution to faster moving staus. Though this does reduce the total number of 
slow-moving staus, the stopping length of the stau depends linearly on its mass. In other 
words, though there are fewer slow-moving staus, they are more easily stopped. We find 
that when both effects are taken into account, the effective number of stopped staus at a 
given distance from the collision point increases as the stau mass decreases. For example, 
comparing the benchmark Maj^jj-, scenario with Apg = 140 GeV to the benchmark Maj^j 
scenario with Apq = 209 GeV, we find the number of stopped staus at a given distance is 
larger in the Majj^j model by a factor of ~ 1.3. 

We now discuss how such stopped staus could in practice be observed. First consider the 
case of staus stuck in the rock surrounding the detector. In this case, a full reconstruction 
of the timing of the stau decay with that of the initial production can likely not be achieved, 
but the observation of the decay product, given by a leptonically decaying tau could still 
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Figure 19: Plot of the number of stopped staus expected from 100 fb^^ of simulated LHC 
data in two F-theory GUT scenarios as a function of distance from the collision point in an 
idealization where surrounding the collision point there is a large ball of either carbon (C) 
or iron (Fe). As we integrate over larger distances, the number of stopped staus necessarily 
increases. 
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be observed. The decay of a tau to a muon can be measured by the muon chamber, and in 
some cases also the inner tracking chamber. 

The primary source of Standard Model background for such muon events is from cosmic 
rays. Such cosmic rays enter from the sky above, and so will correspond to muons passing 
downwards through the detector. By contrast, there are no cosmic rays passing through 
the other side of the Earth, so that contamination from upwards moving cosmic rays will be 
negligible. Hence, we focus on the prospects for detecting upwards moving muons from stau 
decays. The number of such upward moving muons is estimated to be ~ 0.7% of the total 
number of stopped staus 86 . Assuming the number of stopped staus can be estimated by 
the number of staus with /3 < 0.5, we can expect C(50) events in the benchmark Maj^Jjj-, 
scenario. We note that this is an order of magnitude estimate, and should be improved 
with a more realistic simulation. 

Even for upwards moving muons, however, there is still a non-negligible background 
from atmospheric neutrinos which pass through the Earth. When these neutrinos (and 
anti-neutrinos) strike nuclei underneath the detector, they can generate a conversion z/ — )■ /x, 
with roughly equal likelihood for generating /x"*" and Using the results of the MACRO 



experiment 96 , the number of such upward-going muons passing through a 20m x 20m 
detector is expected to be ~ 80 events/year. Using the same selection cuts on upward 
passing muons as in [96], it was found there that the muon energy distribution from this 
background has a broad peak around 50 GeV. The kinematic upper bound on the energy of 
a muon generated by the decay ri — )■ Gi/^i/^/i will have energy < ~ 50 — 100 GeV, 

which is itself on the tail of the distribution of this three-body process. This means that of 
the roughly 80 background muons, less than half actually constitute background for stau 
decays. Further, though this background is not completely negligible, the signal can still be 
statistically significant if the velocity and direction of the background is reliably calculated. 
In addition, the tracking information of a stau event can be helpful in identifying these 
events (86|. This is because the momentum and the direction of the slow-moving staus 
can be recorded during their escape from the detector. This information can be used to 
infer the trajectory along which a slow-moving stau could be stopped. Looking for muons 
coming from a particular direction, this would provide a means to significantly reduce the 
background of such events. 

Staus can also become stuck in the detector. For tau leptons which decay hadronically, 
the optimal location for a stau to be stopped would be in the hadronic calorimeter of 
either the ATLAS or CMS detector. This type of search is similar in spirit to looking 
for stopped gluinos in split supersymmetry models [97]. From figure 19 we see there are 
~ 3000 (resp. 5000) staus stopped in 2 meters (resp. 4 meters) of iron out from the 
collision point for the low A model. This implies that a similar number of staus should be 
stopped in the calorimeter. Let us first consider the case that those staus stopped inside 
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the hadronic calorimeter (HCAL) decay hadronically via the tau decay product of the stau. 
This possibihty is interesting since the tau jet is typically narrow, and therefore can deposit 
all its energy in one or very few calorimeter cells. For a 200 GeV mass stau, the tau jet 
will have energy ~ 100 GeV. In this case, the hadronic calorimeter would light up, but 
no track would be seen in the tracking chamber, and no energy would be deposited in the 
electromagnetic calorimeter (ECAL). This is in contrast to QCD jets which would leave 
tracks as well as have energy deposited in both the ECAL and HCAL. This would seem to 
indicate that QCD background can be reduced significantly. Even so, given the large number 
of QCD events, any leakage to the signal could be dangerous, which must be studied in 
detail separately. Cosmic rays are another potential source of background. However, most 
of the cosmic rays are muons and can be detected in the muon chamber. If they are neutral 
particles, they are likely to be stopped by the outer material of the detector. Stopped 
staus stuck in the HCAL can also decay into muons, although the chances are smaller at 
around 10%. If the muon goes outward from the hadronic calorimeter, the muon can only 
be detected in the muon chamber when there is no activity in other parts of the detector. 
Standard Model objects cannot mimic this type of process, and so hardware fluctuations 
or malfunctions would appear to provide the main source of background. 

Stopped staus can in principle also be used to deduce the scale of supersymmetry break- 



ing [86 93 95 . Recall that the decay rate of the stau NLSP to a gravitino and tau T^^ is 
given in terms of the scale of supersymmetry breaking a/F^ felt by the gravitino and the 
mass of the stau through the relation (see for example [74]): 



rn = (T7^;^) (^^^V-(2xlO-^eV). (40) 



100 GeV/ V VFo J 

Measuring the scale a/Fq would then allow us to also determine the mass of the gravitino. 
Let us note that although here we have distinguished between and the scale of super- 
symmetry breaking a/F felt by the messengers, in the context of F-theory GUTs it is most 
natural to set F = Fq. In F-theory GUTs, F ~ 10^^ GeV^, and the hfetime of the NLSP is 



on the order of one second to an hour 12 



Methods to measure the lifetime of the stopped stau typically rely on being able to 
transport the material in which the stau is stopped to a remote location, where its decay 



can then be carefully studied. For example, in 93 , 95 it was proposed that a "stopper 
detector" be added which could surround the ordinary detector, and in [94j it was proposed 
to use a large water tank as a collector for stopped staus. Staus which become lodged in the 



neighboring cavern rock could in principle also be removed 86 . Both methods work best 
for staus which are stable on the order of at least a day to a week. In the context of F-theory 
GUTs, however, where the stau lifetime is far less than a day, such detection methods do not 
appear as practical. Other methods have been proposed for detecting quasi-stable staus 
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which decay on shorter timescales, as for example in |98|jj For example, in the related 
context of stopped gluino scenarios , it was proposed that if a pair of gluinos become 
stopped in the detector or neighboring cavern, correlating the time separation between the 
decays of these pairs would then provide another means to extract the lifetime, and could 
in principle be adapted to the case of stopped stau scenarios. 

Though a full study is beyond the scope of the present paper, we now speculate on 
potential ways to measure the lifetime of staus in the context of F-theory GUTs. The 
primary challenge in extracting the lifetime of the stau stems from correlating the initial 
production of the stau with its final decay. Although staus typically decay out of time 
with the bunch crossing of the original production event, the information on the stau (for 
example information of tracking and energy deposit in both the HCAL and ECAL) may 
already be recorded in the data if the production event is hard enough (e.g. with hard jets 
or leptons) to be triggered on. In fortuitous circumstances, the direction and velocity of the 
stau could then be reconstructed from an off-line analysis from the energy deposit dE/dx 
in the calorimeters and monitored drift tube data for staus stopped in the muon chamber. 
In such cases, it would then be possible to deduce where such a stau would possibly be 
stopped. For example, if the stau stopped in the middle of the CMS muon chambeiQ the 
tracking apparatus would detect a muon track traversing and then ending in some section of 
the muon chamber, and it would become possible to estimate the position of the candidate 
stopped stau. If at some later time a muon was emitted in the vicinity of the stopped 
stau, and provided there is no other activity in that part of the detector, this would likely 
correspond to a signal for a stopped stau. The time difference between the decay event and 
production event would then provide a measure of the lifetime of the long-lived stau. For 
staus stuck in the HCAL a similar type of correlation analysis may also be available. Given 
the tantalizing prospect of obtaining direct access to the scale of supersymmetry breaking, 
it would seem important to investigate this issue further. 



7 Mass Spectra and Distinguishing Look-Alikes 

Our focus in much of this paper has been on signatures of stau NLSP F-theory GUTs. 
Assuming evidence has been found for such a scenario, it is important to study what 
classes of models reproduce the same observables. This broader question based on the LHC 

^Though not the case studied here, for staus which decay with an average time of 10^'* — 10^'"^ seconds, 
it is possible to exploit the fact that there is a distribution of timescales for the decay of the stau, so that 
some staus will decay inside the detector pO]. For F-theory GUTs, this method is not available because 
the stau has a lifetime of one second to an hour. 

^Compared to ATLAS, the CMS muon tracking system has a far higher density of iron, making it much 
more likely that a stau can become trapped there. We thank members of the CMS group for discussions 
on this point. 
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inverse problem has been studied from different vantagepoints for example in 99 , 100 



One approach to this problem is based on the general footprint method 101 , 102 , which 



directly compares the collider observables of different classes of models to determine possible 
differences. This method has been applied to F-theory GUT models with a bino NLSP 
in 12 to compare particular classes of well-motivated supersymmetry breaking scenarios. 



In this section we consider a simplified version of this problem by studying possible 
degeneracies in the mass spectrum of F-theory GUTs with other models with an MSSM- 
like spectrum. More precisely, we study how mass measurements extracted from stau NLSP 
scenarios can be used to distinguish between F-theory GUTs and minimal gauge mediated 
supersymmetry breaking (mGMSB) models. To this end, we first review how to extract 
details of the sparticle mass spectrum in stau NLSP scenarios which are similar to F- 
theory GUTs. Next, we use the high precision of these mass measurements as a means to 
distinguish between F-theory GUTs and mGMSB models. 



7.1 Mass Reconstruction with Staus 



Following j62j and 63 , we now discuss how quasi-stable stau scenarios allow one to deduce 
much of the sparticle mass spectrum. We can extract the mass of the stau through the 
relation: 

where (5 is measured with the time-of-flight in the precision muon chambers, and p is 



reconstructed through a measurement of pt and rj 88 . Due to the uncertainty in the 



measurement of p and /3, the reconstructed mass distribution will show a broadening around 



the correct mass. From the resolution of p given in equation (37) of section |4| we see that 
the uncertainty in p is proportional to p^ for large stau momentum. This illustrates that 
the mass reconstruction is more accurate at low momenta and velocities. Since the relative 
number of low velocity staus is lower compared with the number of high velocity staus, 
increasing the total number of such events requires a relatively large integrated luminosity. 
It was shown in [62] that the stau mass can be determined up to an uncertainty of ~ 0.02 
GeV using 30 fb~^ of integrated luminosity. 

Reasonable estimates of the stau mass can also be obtained for lower integrated lumi- 
nosity. For example, at ^/s = 14 TeV and ~ 500 pb~^ of simulated LHC data, figure [S 
of section ^ illustrates that for the Maj^Jlj-, benchmark model, the total cross section for 
supersymmetric production of events is ~ 1 pb, indicating there will be about 500 super- 
symmetric events with staus. In this case, we find that the best fit of the stau mass is 
m^** = 175.59 ± 0.47 GeV using a Gaussian fit (see figure 20). In the analysis, we selected 



staus with velocity 0.6 < j3 < 0.91. With 100 pb ^ data, we find only 69 event stau candi- 
dates remain, but the fit of the stau mass is almost as good as the previous case, with the 
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Figure 20: Distribution of the reconstructed stau mass for the benchmark model Maj^jj-, 
defined in subsection 3^ 500 events are generated and staus are selected with velocity 
0.6 < (3 < 0.91. The Gaussian fit gives m^** = 175.59 ± 0.47 GeV, while the actual stau 
mass is 175 GeV. 



error bar increased slightly to 0.63 GeV. 

The presence of quasi-stable staus in the final state provides a convenient means to 
reconstruct the masses of other sparticles further up a decay chain. For example, the mass 
of Xi can be reconstructed by exploiting the fact that over the entire parameter space of 
stau NLSP F-theory GUT scenarios, there is a significant branching fraction to Xi ~^ '^i'^- 
By selecting final states which contain opposite sign staus and two opposite sign r's, we 
can reconstruct the resonance associated to the invariant mass of a pair of rj^r^ to a 
Note that a similar resonance will also occur for xH-, providing a means to reconstruct the 



mass of this sparticle as well. In 62 this type of strategy was employed to work backwards 
up cascade decays to deduce the masses of many sparticles. In particular, the "e model" 
considered in [621 has a similar mass spectrum to that of the Maj^Jlj-, benchmark model 
(though different enough to be distinguishable!), and so rather than perform a complete 
analysis of our case, we shall simply use the same precision estimates based on 30 fb~^ of 
integrated luminosity at ^/s = 14 TeV: 



Am^, = 0.021 GeV, Am^^ = 1.2 GeV, Amj^ = 2.0 GeV 
Am^o = 0.9 GeV, Am^o = 2.0 GeV, 

Xl ' A2 ' 

Arriq^ = 2.8 GeV, Arrig^ = 3.7 GeV, Am^ = 57.7 GeV. (42) 
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Note that for all of the non-colored sparticles, the accuracy is within the range of 2 GeV 
or less. In addition, the mass difference for first and second generation squarks are both 
less then 4 GeV. The reason for the larger uncertainty in the sbottom mass is that events 
from bi and 62 are mixed together since both of them can decay to bxi [62j . The result only 
serves as an estimate of the combination of 61 and 62 masses. 

To distinguish between models which can mimic the spectrum of F-theory GUTs, we 
will also find it convenient to reconstruct the mass of the first and second generation right- 
handed sleptons Ir. In this case, the decay chain of interest depends on the relative masses 
of Xi and Ir. Indeed, we have seen that for A^'io = 1 scenarios at moderate PQ deformation 
mj^ > m^o, whereas at larger PQ deformation for A'^io = 1 and for all values of Apq in 
NiQ = 2 scenarios, we instead have mj^ < m:^o. The latter case has been studied in detail 
in 63 where it was found that with 100 fb~^ of integrated luminosity at y/s = 14 TeV, 



Amr ~ 1 GeV. 

Let us now return to scenarios where mj^ > m^o so that right-handed sleptons domi- 
nantly decay via 

Tr ^ Xil ^ r±r^/. (43) 

Since the actual production of right-handed sleptons is somewhat low in models such as 
the Maj^Jlj-, benchmark model, generating a sufficient number of events requires a large 
integrated luminosity, in part because heavier sparticles rarely decay to Ir, and because the 
direct production cross section <Jj+j- = 1.3 x 10^^ pb is relatively small. To demonstrate 
the reconstruction of the right-handed slepton mass can be done with small enough uncer- 
tainty, we perform an analysis with 100 fb~^ integrated luminosity of simulated data for the 
benchmark model Majj^Jlj-,. First, we selected events with at least one stau candidate, one 
tau jet and no jets. Then the invariant mass Miny{T^T^) of any opposite-sign tau-stau pair 
is calculated]^ Only if there is one invariant mass within 5 GeV around the true mass of x?, 
the event is selected. Finally, we take all possible combinations of leptons in the event with 
the tau-stau pair, and the invariant mass of the lepton-tau-stau system is calculated, giving 



rise to the distribution in figure 21 Fitting the peaks with a Gaussian plus a constant 
background 



Aexp 



1 



B (44) 



^,(M„„(r±r^/)-MP^'^'=)^ 

with A, B, a and M'^'^"'^ fitting parameters, we find that the peak position is 215.8±0.6 GeV. 
This is in good agreement with the actual Ir mass 214.3 GeV. Of course, we should also 
include uncertainties which propagate from the mass measurement of fi and x?- Since these 
uncertainties are typically < 1 GeV, including them does not significantly affect the overall 
uncertainty in the mass measurement of Ir. 
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Here we assume that the tau four momentum can be fuUy constructed. In practice, a cahbration of 



the tau momentum can be done as in 



62 
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Figure 21: Plot of the lepton-tau-stau invariant mass in the benchmark Maj[J|j3 model using 



the event selections discussed below equation (44). The peak at ~ 215 GeV corresponds to 



the resonance from decays of the right-handed selectron and smuon Jlji. The somewhat 
lower and broader peak at ~ 380 GeV corresponds to the resonance from decays of the 
left-handed selectron cl and smuon JIl- 



Let us note that in the above we have quoted results from the literature which typically 
utilize integrated luminosities on the order of 30 — 100 fb~^. We have also seen that 
cruder mass reconstructions can be performed at lower integrated luminosity. It would be 
interesting to study the overall expected mass resolution of such stau NLSP scenarios as 
a function of integrated luminosity. Further detailed information on mass measurements 
would also be readily available with the ILC. 



7.2 Comparison with Minimal Gauge Mediation 

Measuring the masses of the sparticles provides a very precise way to potentially distinguish 
models with similar spectra. Indeed, fitting to an F-theory GUT model leads to a quite spe- 
cific range of possible values for Nio and A. Moreover, because the PQ charge assignments 
are different in the Majorana and Dirac neutrino scenarios (see section [2]), measuring the 
masses of the sleptons and fitting to the rest of the spectrum would allow us to distinguish 
between these two possibilities. This would then provide even more precise information on 
both the size of Apg, as well as the PQ charge assignments of the sparticles. 

In this subsection we consider minimal gauge mediated supersymmetry breaking (mGMSB) 
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models such that the hghtest bosonic and fermionic superpartners respectively given by a 
quasi-stable ri and Xi have the same mass as their F-theory counterparts. We analyze to 
what extent such look-alikes can mimic the rest of an F-theory GUT spectrum. In mini- 
mal gauge mediation, the available parameters which can be varied are the supersymmetry 
breaking scale a/F in the messenger sector, the analogous scale y/F\) which generates the 
mass of the gravitinoj^ the messenger scale Mmessi the effective number of messengers in 
the 5 © 5 A^5, and tan/3. Here we assume that the mass of the messengers are determined 
through the supersymmetry breaking vev of a single GUT singlet. The parameters F and 
Mmess fix the characteristic scale A = F/Mmess- Parameters such as /i and Bfi are deter- 
mined by the condition that proper electroweak symmetry breaking is achieved for a given 
value of tan /3. Let us note that a much broader range of sparticle masses can be achieved 



in more general gauge mediation setups. See for example 103 104 for recent discussion. 

As a first step in identifying potential mGMSB look-alikes, we must require that the 
stau is the NLSP, and that it is quasi-stable on timescales probed by the LHC. Assuming 



that vF ~ equation (40) now implies that vF > 10^ GeV. In particular, this implies 
that the ratio F/M^^,, = AVF < 0.01. Hence, we can reliably use an approximation 
for the mGMSB candidate look-alike in which all soft masses at the messenger scale are 
proportional to A. Later in this subsection we shall consider a potential look-alike where 
we drop this assumption and allow lower messenger scale models, but we will again see that 
this candidate does not effectively mimic F-theory GUTs. 

In quasi-stable mGMSB models with a/F ~ y/F\), the primary deviation between F- 
theory GUTs and minimal gauge mediation models comes from the specific scale of super- 
symmetry breaking, and the effects of the PQ deformation on the soft scalar masses. As a 
first step towards finding potential look-alikes, we focus on the masses of sparticles where 
the effects of the PQ deformation are largely unimportant, in particular the gauginos and 
the first and second generation squarks. Assuming that we have measured enough of the 
gaugino masses, we therefore require 

3NiPA^^^ ~ ATf ^A^'"), (45) 

where the superscript "F" denotes the F-theory model and "m" is for the "mimic" . Next 
consider the masses of the squarks. There is little left-right mixing for the first two gener- 
ations of left and right-handed squarks, and so we focus on this case here. As for all scalar 



^As noted previously, in the context of F-theory GUTs, it is most natural to set F — Fq. 
^"Note that in general, there will be some mixing between the neutral SU{2)l x U{1)y gauginos and the 
neutral Higgsinos, and the charged gauginos and the charged Higgisinos. This can in principle complicate 
the identification of the gaugino mass measurement, though in the present class of models, this is a small 
effect which we shall ignore. 
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partners in niGMSB models, the masses of the squarks at the messenger scale are: 



m;,= —C~,^|Nt^^^"'\ (46) 



a 

^ 47r 



where we have again schematically indicated the dependence on the various fine structure 
constants and quadratic Casimirs of representations through a and Cq. The possible ways 
that an F-theory GUT could reproduce the same squark masses but with a different number 
of messengers and A is by turning on ApQ, and adjusting the messenger scale. As we 
have already seen in section [3| the effects of the PQ deformation are far milder for heavy 
squarks. In particular, the overall shift in the mass of the squarks is on the order of 
^PQ/rrig ~ 0(100 GeV)VO(l TeV) ~ 0(10) GeV, which is a very slight shift. Varying 
the messenger scale Mmess produces logarithmic contribution to the mass-squared, and in 
minimal gauge mediation models, only shifts the squark masses by an effect which is at 
most on the order of O(IOO) GeV. This yields the additional constraint: 



4^C'9V SAtJJ^A^^) = ^Cq^J Nt~^ K^"^^ + O(IOO) GeV. (47) 

Since the squark masses are on the order of a TeV, we conclude that the messenger number 
and A must agree: 

Sivli'^ = Ni'^\ A(^) ~ A^'"). (48) 

Fixing the parameters A^5 and A to agree with the F-theory GUT model allows it to 
mimic the masses of the gauginos and the first and second generation squarks. We now 
discuss whether it is possible to adjust the values of the remaining parameters Mmess and 
tan /3 to exactly match the rest of the spectrum. In the class of models we are considering, 
the effects of tan/? affect the third generation of squarks and sleptons by inducing more 
mixing between left and right-handed sparticles. To isolate the effects of Mmess, we therefore 
focus on the masses of the left and right-handed sleptons cl and cr. Our strategy will be to 
compare the needed Messenger scale for the two soft mass relations in the presence of the 
PQ deformation. We find that at most one of these masses can be matched by adjusting 

-'■''''■mess- 

To arrive at this result, we now study in more precise terms the masses of the left and 
right-handed selectrons as a function of the messenger scale. Letting Ai'^{t) denote the soft 
masses of a minimal gauge mediation model with effective messengers and characteristic 
scale A evaluated at a renormalization group (RG) time t = log(M^g^g/Q^), the condition 
for the selectrons of the F-theory GUT and look-alike to have the same masses is: 

M|(tM)=>l|(t(^)) + g^.A|Q, (49) 
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Table 1: Mass spectra of F-theory benchmark model Maj^J^ and some candidate mGMSB 
look-alike models with A^s = 3 messengers. All masses are in GeV, and we have taken the 
top mass to be 172.4 GeV in our calculation. 



parameter 




mGMSBl 


mGMSB2 


M 

■'■'■'■mess 






Z X lU 


y/F 


2.2 X 10^ 


2.2 X 10^ 


10^ 


tan P 


24.05 


34.7 


24.5 


m~g 


1112 


1113 


1116 




198.6 


199.0 


199.3 


X2 


377.1 


379.4 


378.0 


"^xf 


380.3 


382.3 


381.2 


1106 


1112 


1102 




1059 


1066 


1063 


T^h 


857.6 


866.7 


898.1 




1050 


1047 


1059 




997.2 


982.2 


1014 




1032 


1032 


1046 




383.0 
372.5 


421.7 
412.1 


382.2 
371.6 


ruf^ 


214.3 
175.0 


246.9 
174.8 


204.9 
174.7 




366.1 


400.4 


367.7 




384.0 


422.3 


385.1 


rrih 


114.3 


114.3 


113.8 


rriA 


693.1 


614.2 


623.4 



where t^^'' and t^'^^ respectively correspond to the RG time evolution from the messenger 
scale to the weak scale for the F-theory GUT and the mimic. Here we are implicitly using 
the fact that the Yukawas of the first two generations are small, so that the PQ deformation 
can be treated as an additive shift for / — cr and ej,. Eliminating Apg from this system of 
equations, we therefore find the condition: 



(50) 



QsR MUm) - Mliti"^))- 
In Appendix B we compute the leading order behavior of the righthand side of this equation. 
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the end result of which is: 

qe^ 1 + A^5 9 (a2 



qe^ l + 5A^5/33 220 V"! 



(51) 



In our conventions the weak scale values of the couplings are ai ~ 1/100 and a2 ~ 1/30 
which implies: 

^ ~ 1.5 X , , 52 

gen l + 5iV5/33' ^ ' 

which for = 3 and 6 yields the condition: 

ATg = 3 ^ ^ ~ 4, (53) 
ATg = 6 ^ ^ ~ 6. (54) 

In other words, since the actual ratio of g's is 2 for Majorana scenarios, and 1 for Dirac 
scenarios, a candidate look-alike cannot be found where the masses exactly agree. For 
example, in an F-theory GUT Majorana scenario with A^io = 1, this would lead to a 
mismatch by a factor of 2 in the g's. 

In principle, uncertainties in the mass measurement could still allow a candidate look- 
alike to remain compatible with an F-theory GUT, up to experimental error. In the present 
context, we can view this as reflecting an uncertainty in determining the ratio of the g's. 
Since the mimics of A'^io = 1 Majorana scenarios lead to a smaller deviation from the actual 
ratio of the g's, we focus on this case. Rescaling q by the factor of 2 discussed below 



equation (54) translates to a shift in the mass squared of order 2 x Apg. We now ask 
whether this type of shift could be explained by an uncertainty in the mass measurement. 
Using the relation: 

^actual ^ ^observed ^ ^55^ 

and squaring, we obtain: 

^^acW^2 _ ^^observedY ± 2 X m X Am. (56) 

In other words F-theory GUTs which can be mimicked have A pq which satisfies the upper 
bound: 

2 X A%Q < 2 X m X Am. (57) 

In typical F-theory GUTs, the masses of the right-handed sleptons satisfy m ~ 300 GeV. 
Combined with the uncertainty in the mass measurement Am < 2 GeV reviewed in the 
previous subsection, we find: 

ApQ < 30 GeV, (58) 
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Table 2: Mass spectra of F-theory benchmark model Majjjf and some candidate mGMSB 
look-alike models with = 3 messengers. All masses are in GeV, and we have taken the 
top mass to be 172.4 GeV in our calculation. 



parameter 


MajS 


mGMSB3 


mGMSB4 




10 


10* 


3 X 10^ 


\/F 


2.2 X 10^ 


2.2 X 10^ 


1.22 X 10^ 


tan /3 


24.9 


40.2 


40.7 


m-g 


1110 


1118 


1128 




198.3 


199.7 


200.0 


373.8 


376.1 


350.9 




377.5 


378.7 


350.3 


1098 


1098 


1091 




1050 


1062 


1062 




846.1 


916.3 


968.4 




1041 


1049 


1062 




979.8 


986.3 


1008 




1019 


1037 


1052 




329.7 
317.4 


363.7 
352.9 


330.0 
318.4 




166.6 
104.6 


188.9 
104.5 


163.3 
104.7 


mo. 


309.5 


345.2 


314.5 




331.5 


373.2 


340.5 




114.3 


113.3 


112.8 


rriA 


701.8 


497.5 


395.9 



which should be viewed as an order of magnitude estimate. In other words, the high 
precision of the mass measurement translates into a tight upper bound on the size of the 
PQ deformation which can be mimicked. 

To further support this general analysis, we now consider in greater detail potential 
look-alikes for the benchmark model Maj^j^. As a first example, consider mGMSB models 
with Mmess = 10^^ GeV, A = 50 TeV, and = 3. The only remaining parameter of an 
mGMSB model which we can vary is tan f3, which changes the amount of left-right mixing 
for staus, allowing the lightest stau to become degenerate in mass with the F-theory GUT 
value. By increasing tan/3 to 34.7, one can achieve a stau NLSP with mass very close to 
that of the F-theory GUT. The mass spectrum of this model (mGMSBl) is shown in Table 
[Tj Though the gaugino and squark masses are within a few GeV, the slepton masses in 
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mGMSBl deviate from the Maj[J|j-, values by as much as ~ 30 GeV. Returning to equation 



(42), we conclude that this difference is large enough to distinguish these two models with 
sufficient LHC data. 

Next consider variations of the messenger scale Mmess- Table [T] displays the mass spec- 
trum of an mGMSB model (mGMSB2) where the model parameters are the same as the 
benchmark Maj^jj-, except Mmess = 2 x 10^ GeV and tan/3 = 24.5. One can see that in 
addition to matching the gauginos and first and second generation squark masses, the left- 
handed slepton masses are also very close to the Maj^jj-, model. The notable differences 
are the right-handed slepton mass and the third-generation squarks, as well as the Higgs 
masses. The difference in the third-generation squark masses is mainly due to the lower 
messenger scale which shortens the contribution from renormalization group flow effects. 

For F-theory models with larger Apg, the deformation effects are larger. As we have 
seen, this makes it more difficult to find an mGMSB look-alike. As an example we consider 
the benchmark model Maj[|j where A'^io = 1, A = 50 TeV and Apg = 209 GeV. Even when 
the messenger scale is as low as ~ 10*^ GeV in the model mGMSB3, table |2] shows that 
the slepton masses typically differ by as much as ~ 35 GeV, which can be resolved with 
the mass measurement procedure detailed above. For even lower messenger scale, the stau 
NLSP in mGMSB model would decay inside the detector, so long as we assume Fq = F. 

Even when we allow the possibility that Fq ^ F so that the stau decays outside the 
detector, we find that representative low scale messenger models do not correctly mimic the 
spectrum of the Maj^]^ scenario. An example of this type is shown as the model mGMSB4 
in table |2] Let us note that when F/M^g^^ ~ 0(1), the theoretical analysis presented 
earlier in this subsection and in Appendix C is no longer accurate. Table [2] shows that even 
though the slepton masses can be quite degenerate with the above mentioned F-theory 
model, the masses of xf would be ~ 20 GeV lighter than the corresponding masses 

in the F-theory model. The reason for this is mainly due to the lowered messenger scale 
where the running of m|^^ is reduced and therefore become less negative. From the relation 

l/^P ~ -"^k - ^"^1' (59) 

we can see the n term decreases in this situation, and therefore the higgsino masses be- 
come closer to the gaugino masses. The increased mixing between gauginos and higgsinos 
decreases the physical masses of ^2 xf- This mismatch in the neutralino and chargino 
masses is by itself enough to resolve the potential degeneracy. 

Our numerical analysis has mainly been based on two benchmark models with A^io = 1. 
The basic idea is similar for other models, as well as for A^io = 2 scenarios. Though beyond 
the scope of the present paper, it would be interesting to perform a more complete scan 
over other well-motivated look-alikes. 
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8 Conclusions 



In this paper we have seen that over a broad range of parameter space in minimal F- 
theory GUTs, a quasi-stable stau is the NLSP. The consequences for collider physics are 
quite striking, since in stau NLSP scenarios, this particle will appear as a massive quasi- 
stable charged particle at the LHC. We have seen that with the center-of-mass energies and 
integrated luminosity expected at the LHC, signatures of this scenario can be discovered. 
We have also seen how the corresponding number of stau events correlates with the number 
of leptonic final states as a function of F-theory GUT parameters. Focussing on inclusive one 
and two stau events, it is possible to constrain the characteristic scale of gauge mediation 
A, and using more exclusive signatures such as Drell-Yan production of staus, it is possible 
to constrain the size of the stringy PQ deformation parameter present in F-theory GUTs. 
Slow-moving staus can also become stuck in or near the detector, and we have analyzed in 
an idealized situation the number of such events to expect. Finally, we have used the high 
precision expected in mass measurements in stau scenarios as a way to compare the spectra 
of F-theory GUTs with minimal gauge mediation models. 

In the remainder of this section, we briefly speculate on some possible avenues of further 
investigation. In this work we have focussed on the potential for discovering quasi-stable 
stau NLSP F-theory GUT scenarios. Quasi-stable stau NLSPs are also a feature of other 
models with an MSSM-like spectrum, and we have seen that there are some measurable 
differences between minimal gauge mediation models and F-theory GUTs. Broadening the 
range of possible models both to more general gauge mediation models, and more general 
scenarios with a stau NLSP, it would be interesting to further study the level at which 
F-theory GUTs can be distinguished from other well-motivated possibilities. 

We have seen that in the context of F-theory GUTs, staus can be become stuck in or near 
the detector. This holds out the exciting prospect of measuring the scale of supersymmetry 
breaking. It would be interesting to further investigate with more realistic simulations how 
to use the ATLAS and CMS detectors to detect stopped staus, and measure their lifetimes. 
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Appendix A: Simulation Details 



We simulate SUSY events using PYTHIA 83 . The analysis is done based on the parton-level 
objects 

• lepton: pt > 5 GeV and \ri\ < 2.5 and the transverse energy deposited in a cone 
AR= 0.4 around the lepton is less than 5 GeV. Note: staus which have large velocity 
P > 0.91 are also treated as muons. 

• tau jet: pr > 20 GeV and \ri\ < 3.0, with efficiency = 0.5. 

leptonic decaying taus are identified as leptons based on the lepton cuts. 

• jet: Pt > 20 GeV and \r]\ < 3.0. 

We shall often augment these basic cuts by additional analysis cuts. 

We include the detector resolution of the momentum and energy for muon, lepton and 



jet, which are parameterized as follows in the ATLAS detector 88 
• For muon the momentum resolution is: 



Pt TeV ^/sm6 
Here we take a = 36% and b = 1.3%, 
For electron the energy resolution is: 

aE a 



(60) 



(61) 



Here we take a = 10% and b = 0.4%, 

For jet, the energy resolution is similar to electron, but with a = 80% and b = 15%. 



Appendix B: Selectron Mass Comparison 



In this Appendix we compare in more precise terms the masses of the left and right- 
handed selectrons in F-theory GUTs and minimal gauge mediated supersymmetry breaking. 
In minimal gauge mediation, the first and second generation soft masses squared at a 
renormalization group (RG) time t = log(M^g^^/(5^) are given by (see for example |74]): 



ar(t) 
47r 
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An 
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- 1 



(62) 
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where we have dropped running effects due to the Yukawa couphngs, which are small for 
the first two generations. In the above, r is an index labelling the three g auge groups, is 
the quadratic Casimir associated to the representation of the scalar, ki = b/^, k2 = = 1 
so that ai{t)ki unify at the GUT scale, and in this normalization br — 11, 1, —3 denote the 
beta functions for r — 1,2, 3, respectively. Further, the renormalized fine structure constant 
satisfies: 



0:^(0) = ar{t) 



4:71 



brt 



(63) 



Here, t = corresponds to the Messenger scale, and we shall be interested in evaluating 
both the F-theory GUT and candidate look-alike at the weak scale. Because of this, we can 
replace each ar{t) by its weak scale value, which we denote by q;^. In this case, all of the 
difference from RG flow effects is captured by the terms in brackets containing arbrt/A7r. 
Letting t^^^ and t'^'^^ denote the respective values of the RG time for the F-theory GUT 
and the candidate look-alike, it follows that the right-handed selectrons and sleptons of a 
candidate mimic must satisfy the relations: 
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(64) 

For the right-handed selectrons, the sum on r truncates to the term r — 1, and for the left- 
handed selectrons, the r — 2 term dominates. Taking the ratio of these two contributions, 
we find: 



Qe, 1 + 57V5/33 A;iCf 
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(65) 

(66) 
(67) 



where in the second to last equality we have kept only the leading logarithms. Note that 
this final expression is independent of the logarithmic dependence on the messenger scale, 
as well as A. 



59 



Appendix C: Miscellaneous Plots 

In this Appendix we collect additional plots of potential interest. 

Mass (GeV) 
250 




Figure 22: Plot of the four lightest sparticle masses as a function of messenger scale Mmess- 
Here A'^io = 2, A = 28 TeV, and we have displayed the values for zero PQ deformation 
(upper lines) and Apg = 150 GeV (lower lines). Note that the bino mass remains constant, 
both as a function of messenger scale and of PQ deformation. 
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GeV 

g 




Figure 23: Spectrum of an F-theory GUT Majorana neutrino scenario with A^io = 2, A = 28 
TeV for minimal (red left columns) and maximal (blue right columns) PQ deformation. As 
opposed to the case of A'^io = 1 models, for larger numbers of messengers the stau is always 
the NLSP for F-theory GUTs. 
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Mass (GeV) 




Figure 24: Mass spectrum of the sleptons, lightest chargino and two hghtest neutrahnos 
with NiQ = 2, A = 28 TeV as a function of Apg. The dashed grey hne at 100 GeV denotes 
the present experimental bound on the mass of a quasi-stable stau. 
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0.35 




Figure 25: Dominant branching fractions of Xi ^ function of Apg with A'^io = 2 and 
A = 28 TeV. As opposed to the case A^'io = 1, the overall Apg dependence is smaller. In 
addition, the decay channel to sneutrinos is available even for Apg = 0. 
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Figure 26: Dominant branching fractions of X2 a function of Apq with A'^io = 2 and 
A = 28 TeV. As opposed to the case A^'io = 1, the overall Apg dependence is smaller. In 
addition, the decay channel to sneutrinos is available even for Apg = 0. 
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Figure 27: Contour plot of the inclusive > Ifi signal in parameter space A and Apg for 
^/s — 14 TeV and A^io = 2 messengers. The three contours correspond to the value of the 
cross section at 0.5 pb, 0.05 pb, and 0.01 pb. To detect five inclusive one or two stau events, 
this amounts to an integrated luminosity of respectively 10 pb~^, 100 pb~^ and 0.5 fb~^. 
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V7 (TeV) 

Figure 28: Plot of the leading order cross sections of various production channels as func- 

(2) 

tions of center-of-mass energy ^/s for the benchmark Majjjf scenario. 
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